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PREFACE 


This  project  is  the  responsibility  of  the  Chemistry  Program  of 
the  Physical  Sciences  Division  of  Stanford  Research  InBtltutei  In 
cooperation  with  the  Physics  and  Chemistry  of  Fluids  Program.  Project 
organization  and  principal  contributors  to  the  technical  work  are: 


Program  Manager:  M.  E.  Hill 

Low-Velocity  Detonation:  R.  W.  Woolfolk 

H.  R.  Bredfeldt 
D.  tegg 

Physics  and  Chemistry  of 

Detonation:  R.  Shaw 

M.  Cowperthwaite 

Kinetics  and  Mechanisms 

of  Thermal  Decomposition:  D.  S.  Ross 

D.  Tegg 
M.  E.  Hill 


SUMMARY 


Under  the  sponsorship  of  the  Office  of  Naval  Research,  Stanford 
Research  Institute  is  studying  the  following  fundamental  sensitivity 
properties  of  difluoramino  compounds: 

1.  The  low-velocity  detonation  (LVD)  wave  characteristics 
of  the  liquid  phase, 

2.  The  relation  of  shock  sensitivity  and  failure  diameter 
to  the  flow  and  the  chemical  reaction  rate  behind  the 
shock  front. 

3.  The  mechanism  and  kinetics  of  thermal  decomposition  of 
the  compounds  in  the  gas  phase. 

The  compounds  chosen  for  this  study  are  as  simple  as  possible 
structurally,  but  at  the  same  time  they  have  physical  properties  that 
permit  their  use  in  each  phase  of  the  study.  The  model  compounds  are 
the  bis  difluoramino  and  trJ.s  difluoramino  Isomers  of  the  propane 
series  primarily  1,2-,  1,3-,  1,1-,  and  2, 2-bis(difluoraalno)propane  } 
(1,2-,  1,3-,  1,1-,  and  2,2-DP);  1, 2-bis(dlfluoramino)-2-methyl  propane 
(IBA);  2, 3-bis (difluoramino) butane,  (2,3-DB)- in  meso  and  dl  forms; 

1,2, 2-trls(difluoramino)propane  (1,2,2-H>);  and  2,3,3-tris(difluoramino)- 
butane  (2,3,3-TB).  These  compounds  are  being  used  so  that  the  results 
from  each  experimental  technique  can  be  interrelated  without  having 
to  account  for  effects  caused  by  different  functional  groups,  amount 
of  carbon  content,  and  chain  branching.  Through  the  use  of  isomers, 
structural  effects  on  the  results  from  a  particular  technique  can  be 
determined  for  example,  the  differences  in  the  effect  of  vicinal  and 
geminate  difluoramino  groups  on  the  ease  of  initiation  and  on  the 
environmental  stability  of  compounds  containing  them.  'A  similar 
approach  has  begun  with  isomeric  nltro  compounds. 

Low-Velocity  Detonation 

Hie  shock  wave  interactions  that  can  lead  to  the  initiation  of 
low-velocity  detonations  (LVD)  are  examined  with  particular  emphasis 


v 


on  cavitation.  The  two  forms  of  cavitation,  gas  inclusion  and  pure 
vaporization,  are  discussed.  The  latter's  relationship  to  the  shock¬ 
wave  interactions  in  the  LVD  gap  sensitivity  tests  is  examined.  We 
believe  that  pure  vaporization  cavitation,  which  cannot  be  removed 
by  degassing,  is  present  under  gap  tjst  conditions.  Experimental  studies 
in  the  gap  test  arrangement  show  the  present  of  cavitation  in  qualitative 
agreement  with  theoretical  prediction. 

The  "universal  Hugoniot"  for  liquids  has  been  reevaluated  using  a 
temperature  correction  for  the  sonic  velocity,  c^.  The  new  equation  is 
improved  slightly  over  the  one  formerly  used.  This  concluded  the 
effort  in  this  area. 

A  trace  dye  method  for  following  the  flow  field  behind  a  shock 
wave  was  studied.  We  concluded  that  this  method  was  not  suitable  at 
present  for  our  studies  and  work  on  it  will  be  discontinued. 

Experiments  on  the  effect  of  confinement  wall  thickness  on  the 
LVD  gap  sensitivity  are  discussed  and  preliminary  charge  diameter  studies 
are  reported.  A  summary  of  the  effects  of  confinement  on  the  LVD 
gap  sensitivity  is  presented. 

Some  of  the  physical  and  chemical  properties  of  liquid  explosives 
are  correlated  with  the  LVD  gap  sensitivity.  This  correlation  is 
remarkably  good  considering  its  simplicity.  The  evidence  that  LVD 
initiation  may  be  a  gas  phase  phenomenon  is  discussed.  This  further 
supports  the  contention  that  LVD  may  be  related  to  the  formation  of 
cavities  in  the  liquid  explosive. 

Physics  and  Chemistry  of  Detonation 

A  new  method  for  calculating  shock  temperatures  in  liquid  explosives 
has  been  developed  from  the  Walsh-Christlan  method  for  metals.  The  Walsh- 
Chri.*tian  method  assumes  that  Cv  and  (dp/3T)y  are  constant.  The  basis  of 
the  pr  lent  method  is  a  variational  analysis  which  shows  that  shock  tempera¬ 
ture  more  sensitive  to  variations  in  Cy  than  to  variations  in  (dp/dT>v. 
Thus  Cv  *  1  assumed  to  be  a  function  of  temperature  and  Op/dT>v  is  assumed 


vi 


to  be  constant.  It  is  considered  likely  that  Cy  will  increase  with 
temperature  due  to  an  increase  in  the  vibrational  heat  capacity.  The 
simplest  assumption  is  that  C  of  the  liquid  will  increase  as  C  of 

V  V 

the  ideal  gas.  This  is  the  basis  of  the  Cy(T)  model. 

Shock  temperatures  of  carbon  tetrachloride  and  nltromethane  have 
been  calculated  using  both  constant  Cy  and  Cv(T>.  Shock  temperatures 
for  carbon  tetrachloride  calculated  using  the  CMT)  model  are  in  better 
agreement  with  experimental  brightness  temperatures  than  those  calculated 
using  the  Walsh-Christian  constant  Cy  model.  Shock  temperatures  for 
nltromethane  are  compared  with  those  of  other  workers. 

Shock  temperatures  in  water  were  calculated  using  constant  C  and 

y 

(3p/aT)y.  The  temperatures  are  about  15£  less  than  those  calculated 
by  Rice  and  Walsh  using  a  model  based  on  constant  C  .  This  difference 
is  not  regarded  as  significant  compared  with  the  uncertainties  in  the 
variation  of  Cy  and  with  temperature  and  pressure. 


The  VLPP  studies  of  all  the  NF  compounds  we  have  had  at  hand  have 
been  completed,  and  the  accumulated  data  suggest  a  model  of  unimolecular 
decomposition  in  which  the  NF,  compound 


is  destabilized  by  both  o-NFa  substitution  due  to  a  decrease  in  E  ,  and  by 
B ~alkyl  or 3 -nf2  substitution  due  to  an  increase  in  the  A  factor.  Since 
both  factors  are  operative  in  the  case  of  2,2-DB,  it  is  thermally  the 
least  stable  NF  material  we  have  studied. 

A  detailed  reaction  coordinate  diagram  may  be  constructed  for  the 
system  FaN-C-C-NFa  j!  C=C  +  N,F4  on  the  basis  of  both  our  work  and  that 
by  Trotman-Dlckenson.  It  is  possible,  on  the  basis  of  the  diagram,  to 
explain  the  side-products  observed  in  the  VU>P  of  1,2-DP. 

The  VLPP  of  1-nit rcpropane  has  been  studied.  The  predicts  are 
propylene  and  presumably  HONO,  and  the  T,/a  value  is  665°C. 
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I  INTRODUCTION 


Under  the  sponsorship  of  the  Office  of  Naval  Research,  Stanford 
Research  Institute  is  studying  fundamental  sensitivity  properties  of 
liquid  high-energy  materials  in  order  to  obtain  basic  knowledge  on  the 
ease  of  Initiation  and  the  characteristics  of  propagation,  especially 
in  relationship  to  chemical  structure.  In  this  work  we  have  been 
primarily  concerned  with:  (a)  detonatability  of  the  liquid  phase, 
i.e. ,  whether  it  will  support  a  detonation  wave,  and  the  failure  diameter 
for  detonation;  (b)  the  necessary  conditions  for  initiating  detonation 
and  whether  ease  of  initiation  is  related  to  structure;  and  (c)  the  mode 
of  decomposition  of  the  chemical  structure  and  the  relationship  of  the 
initial  decomposition  steps  with  the  phenomena  of  (a)  and  (b). 

In  previous  work  on  this  contract,1  these  relationships  have  been 
studied  with  isomeric  difluoramino  propanes  and  butaneB  and  we  heve 
now  begun  study  on  nitro  compounds  in  the  analogous  isomeric  propane 
series.  Thus  a  direct  comparison  between  the  two  high-energy  groupings 
which  are  of  most  interest  at  present  cat,  be  made,  and  basic  knowledge 
can  be  obtained  on  the  initiation  and  propagation  properties  of  model 
liquid  nitroaliphatic  plasticizers. 

This  program  consists  of  three  primary  tasks: 

A-  The  detonation  wave  characteristics  of  the  liouid  phase 

Difluoramino  compounds  are  detonatable  at  low  failure  diameters 
and  possess  a  high-velocity  detonation  rate  of  about  6  mm/^sec.  However, 
they  also  have  a  more  easily  initiated  low-velocity  detonation  region. 

The  most  recent  work  has  emphasized  study  of  low-velocity  detonation 
(LVD)  particularly  aimed  at  characterizing  this  unusual  phenomenon 
(Section  II).  Specifically,  the  objectives  are  to  determine  if  the 
phenomenon  is  general  under  conditions  of  larger  diameter  and  differing 
materials  of  confinement,  to  determine  how  easily  LVD  initiation  occurs, 
and  to  relate  the  LVD  to  chemical  decomposition  mechanisms. 
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B.  Phvslca  and  chemistry  of  detonation 

The  ultimate  objective  of  studies  of  the  physics  and  chemistry  of 
detonation  (Section  III)  is  to  correlate  transient  detonation  phenomena, 
such  as  shock  initiation  and  failure  behavior,  with  the  mode  of  decom- 
position.  This  objective  involves:  (a)  measurement  of  failure  diameters; 
(b)  study  of  events  in  the  liquids  as  shocks  of  various  magnitudes 
enter;  and  (c)  measurements  of  reaction  times  at  pressures  and  tempera¬ 
tures  comparable  to  those  encountered  in  initiating  shocks.  During 
the  past  six  months,  reaction  time  measurements  made  in  shocked  difluoramlno 
compounds  and  calculation  of  the  shock  temperature  have  received  con¬ 
siderable  emphasis;  these  measurements  are  pertinent  to  failure  diameter 
theory,  dark  waves,  and  high-temperature  chemical  kinetics  in  the  liquid 
phase. 

C.  The  mechanism  and  kinetics  of  thermal  decomposition  in  the  gas  phase 

Section  IV  of  this  report  comprises  a  study  of  the  decomposition 
behavior  of  NF  compounds  and  describes  a  relationship  between  detonation 
phenomena  and  chemical  mechanisms.  Most  recently  emphasis  has  been 
given  to  obtaining  knowledge  of  the  first  few  steps  of  decomposition 
under  conditions  of  very-low-pressure  pyrolysis  (VU*P)  and  correlating 
these  with  the  low-velocity  and  high-velocity  detonation  phenomena.  In 
addition,  the  VLPP  study  provides  knowledge  of  the  initial  bond-breaking 
reaction  and  bond  energies  of  simple  difluoramlno  groupings.  The  study 
is  being  extended  to  an  isomeric  series  of  dinit ropropane  compounds; 
during  this  period  some  initial  data  on  mononitropropanes  were  obtained. 

The  Appendix  to  this  report  contains  a  list  of  papers,  which 
document  various  phases  of  this  program  and  are  in  the  process  of 
publication  in  the  technical  literature. 
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II  LOW-VELOCITY  DETONATION 


(R.  W.  Woolfolk  and  H.  R.  Bredfeldt) 


A.  Introduction 

Wt  have  been  studying  low-vslocity  detonation  (LVD)  In  the 
difluoramino  compounds.  Our  efforts  have  concentrated  on  the  sensitivity 
of  these  materials  to  LVD  initiation.  The  gap  sensitivity  test  has 
been  the  primary  means  of  determining  the  important  LVD  initiation 
parameters.  Using  the  difluoramino  propanes  and  butanes  and  other 
appropriate  model  compounds  we  have  determined  the  effect  of  confinement 
geometry,  sonic  velocity,  wall  thickness,  and  charge  diameter  on  the 
apparent  LVD  gap  sensitivity.  We  are  now  beginning  to  investigate  the 
importance  of  physical  and  chemical  properties  of  materials  on  LVD 
initiation.  Our  theoretical  and  experimental  studies  have  shown  that 
precursor  cavitation  does  exist  under  gap  test  conditions.  We  have 
found  agreement  between  what  Is  predicted  theoretically  and  what  is 
found  experimentally.  The  presence  of  precursor  cavitation  has  previously 
been  reported  by  the  workers  at  the  Bureau  of  Mines.  1 

Believing  that  cavitation  may  play  a  major  role  in  LVD  initiation 
and  propagation,  we  derived  a  correlation  between  the  heat  of  reaction 
<AHr),  the  heat  of  vaporization  (AlO,  and  the  LVD  gap  sensitivity. 

This  correlation  was  remarkably  good  considering  its  simple  nature. 

The  "universal  Hugoniot"  for  liquids  was  improved  by  the  addition 
of  new  data  and  the  use  of  sonic  velocities  corrected  to  the  experimental 
temperature.  Because  we  feel  further  refinements  are  not  Justified 
by  the  imprecise  nature  of  the  data,  this  completes  our  efforts  in  this 
area. 

We  have  shown  that  the  use  of  2-(2’  ,4 ' -dinit robenzyl) -pyridine 
(DNBP)  to  follow  shock  wave  flow  fields  is  not  practical  under  our  test 
conditions  and  further  work  in  this  area  will  not  be  continued. 
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B.  Theoretical  Studies 


1.  Universal  Hugonlot 

Previous  reports9’1  have  described  an  empirical  method  for  estimating 
the  Hugonlot  of  liquids.  This  equation  was  developed  using  the  data  of 
Rice  and  Walsh4  with  later  addition  of  some  data  from  Cook  and  Rogers.1 
A  recent  Russian  publication9  suggested  the  following  representation 

D  =  1. 2(c  )  +  1.7  u  (1) 

sop 

where  is  the  shock  velocity.  cq  is  the  sonic  velocity,  and  u^  the 
particle  velocity.  Our  representation  (Eq.  2)  differed  from  Eq.  (1)  in 
u  ’<ght,  but  what  appeared  to  be  a  statistically  significant,  manner. 

U  =  1.3,  (c  )  +  1.6.u  (2) 

8  O  P 

Since  both  equations  had  used  some  of  the  same  data  we  did  not 
understand  the  difference.  Therefore  we  undertook  a  review  of  the  data. 
New  results  by  Dick7  of  LASL  were  added  and  the  data  of  Cook  and  Rogers 
were  removed  because  they  did  not  include  their  initial  experimental 
temperature.  We  normalized  the  data  with  cq  corrected  to  the  experimental 
temperature.  Only  those  data  below  the  transition  point  reported  by 
Dick  were  used.  The  temperature  correction  for  cq  was  obtained  from 
the  literature  or  estimated  by  the  methods  outlined  in  Reid  and  Sherwood.9 
The  corrected  data  appear  in  Table  1. 

A  least  squares  fit  of  this  data  (Fig.  1)  is  represented  by  Eq.  (3) 

U  =  1.3.  (c  )  +  1.6.  u  (3) 

s  o  p 

The  "universal  Hugonlot"  in  the  p-v  plane  appears  in  Fig.  2.  The 
curve  through  the  data  was  obtained  using  the  Hugonlot  relationship  and 
Eq.  (3).  The  scatter  in  the  data  was  reduced  by  applying  the  temperature 
correction,  but  some  scatter  remains,  reflecting  the  difficulties  in 
obtaining  very  accurate  experimental  Hugonlot  data  on  liquids.  The 
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Table  1 


SHOCK  PROPERTIES  OP  LIQUIDS 


Compound 

T 

°C 

Co 

mm/|isec 

Po 

g/cc 

P 

kbar 

> 

nun/iisec 

UP 

mn/tiaec 

o 

> 

Benzene 

12 

1.379 

.877 

147 

6.00 

2.75 

.541 

14 

1.367 

.885 

75.3 

4.79 

1.78 

!  .630 

14 

1.367 

.885 

47.1 

4.05 

1.31 

.676 

a16 

1.356 

.883 

52.4 

4.10 

1.45 

.647 

18 

1.345 

.881 

52.1 

4.09 

1.45 

.646 

19 

1.340 

.890 

110 

5.52 

2.26 

.591 

19 

1.340 

.880 

41.3 

3.85 

1.22 

.684 

20 

1.335 

.879 

29.8 

3.47 

.98 

.719 

22 

1.325 

.877 

47.6 

4.05 

1.34 

.669 

22 

1.325 

.877 

15.0 

2.78 

.61 

.776 

24 

1.313 

.875 

106 

5.64 

.21 

.596 

24 

1.313 

.875 

29.1 

3.44 

.97 

.719 

27 

1.298 

.871 

75.4 

4.77 

1.81 

.620 

28 

1.293 

.870 

99.5 

5.28 

2.17 

.589 

28 

1.293 

.870 

88.7 

5.00 

2.04 

.591 

28 

1.293 

.870 

25.8 

3.31 

.90 

.730 

28 

1.293 

.870 

18.9 

2.96 

.73 

.753 

29 

1.288 

.869 

67.4 

4.52 

1.72 

.620 

29 

1.288 

.869 

62.3 

4.38 

1.63 

.627 

29 

1.288 

.869 

17.4 

2.  72 

.73 

.730 

30 

1.283 

.868 

124 

5.71 

2.50 

.562 

a32 

1.273 

.866 

121 

5.66 

2.47 

.564 

Carbon- 

dlsulflde 

a14 

1.175 

1.272 

55.8 

3.47 

1.26 

.636 

17 

1.165 

1.267 

58.5 

3.37 

1.415 

.580 

18 

1.161 

1.266 

61.9 

3.51 

1.39 

.603 

19 

1.159 

1.264 

50.0 

3.39 

1.17 

.655 

20 

1.155 

1.263 

36.5 

3.09 

.94 

.697 

22 

1.149 

1.260 

56.4 

3.47 

1.29 

.628 
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Table  1  (Cont’d) 


Compound 

H 

Cyo 

mm/^sec 

Po 

g/ec 

P 

kbar 

us 

mm/i^sec 

mm  msec 

v/vo 

22 

1.149 

1.260 

18.5 

2.47 

.59 

.761 

24 

1.143 

1.257 

35  7 

3.06 

.93 

.696 

28 

1.130 

1,251 

31.6 

2.94 

.86 

.707 

28 

1.130 

22.9 

2.59 

.71 

.727 

29 

1.126 

21.3 

2.41 

.71 

.705 

32 

1.116 

1.245 

53.6 

3.43 

1.26 

.634 

Carbon- 

14 

.964 

1.606 

106 

4.08 

1.61 

.604 

chloride 

14 

.964 

1.606 

66.7 

3.46 

1.20 

.652 

18 

.956 

1.598 

74.0 

3.50 

1.33 

.621 

19 

.948 

1.596 

153 

4.71 

2.04 

.566 

19 

.948 

1.596 

58.7 

3.28 

1.12 

.658 

20 

.945 

1.594 

42.3 

2.95 

.90 

.694 

a22 

.943 

1.591 

73.9 

3.51 

1.325 

.623 

22 

.943 

1.591 

67.4 

3.44 

1.23 

.643 

22 

.943 

1.591 

21.3 

2.32 

.58 

.752 

24 

.936 

1.586 

148 

4.66 

2.00 

.572 

24 

.936 

1.  586 

41.3 

2.91 

.90 

.692 

27 

.925 

1.581 

106 

4.07 

1.65 

.596 

28 

.921 

1.579 

140 

4.52 

1.96 

.566 

28 

.921 

1.579 

125 

4.27 

1.86 

.565 

28 

.921 

1.579 

36.4 

2.79 

.82 

,702 

28 

.921 

1.579 

26.6 

2.47 

.69 

.722 

29 

.919 

1.577 

95.1 

3.86 

1.57 

.594 

29 

.919 

1.577 

88.1 

3.  74 

1.49 

.601 

32 

.910 

1.571 

24.7 

2.27 

1.69 

.697 

32 

.910 

1.571 

62.8 

3.32 

1.20 

.637 

Acetone8 

26 

1.77 

.7835 

105.8 

5.37 

2.51 

.533 

30 

1.157 

.7793 

46.4 

3.97 

1.49 

.624 

6 


U/c 


Table  1  (Cont ’d) 


Compound 

T 

°C 

Co 

mm  sec 

Po 

e/co 

P 

kbar 

J. 

mm/jiiec 

am/fitec 

v/vo 

Ethyl 

ethera 

21 

.996 

.7122 

QS 

3.88 

1.52 

.610 

32 

.936 

.6995 

E 

5.40 

2.55 

.528 

Ethanol a 

21 

1.158 

.7893 

4.03 

1.487 

.632 

26 

1.140 

.7844 

5.63 

2.50 

.556 

Methanol® 

15 

1.139 

.7960 

46.6 

3.93 

1.483 

.625 

®24 

1.108 

.7877 

109.5 

5.51 

2.525 

.542 

Hexane® 

19 

1.086 

.6599 

41.5 

4.02 

1.517 

.623 

®32 

1.031 

.6487 

95.7 

5.54 

2.59 

.533 

aData  from  Rice  and  Walsh  (Ref.  4). 
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FIGURE  2  NORMALIZED  p-v/vQ  PLOT 
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ugoniots  for  the  difluoroamino  propane  and  butane  that  appeared  in  the 
ast  report  were  calculated  using  Eq.  (2).  The  difference  between 
Eqs.  (2)  and  (3)  is  not  large  enough  to  warrant  a  recalculation  and  will 
have  very  little  effect  on  the  shock  temperature  calculation  pre.ented 
by  Shaw.’  (See  also  Section  111. ,  We  feel  that  the  present  representation 
for  a  universal  Hugoniot"  (Eq.  3)  i.  adequate  for  estimation  purpose, 
and  we  will  not  devote  further  time  to  it.  development. 

2>  ~ — R°le  of  Citation  in  Liquid  Sensitivity 


me  term  "cavitation"  refers  to  the  formation  of  liquid-vapor  or 
liquid-gas  interfaces,  either  at  a  liquid-solid  boundary  or  in  the 
interior  of  a  liquid.  In  the  former  case,  the  resulting  cavities  are 
often  relatively  large,  as  that  found  behind  a  dense  object  dropped  into 
a  pool  of  water  or  around  a  ship-,  propeller.  On  the  other  hand,  cavities 
formed  in  the  interior  of  a  liquid  generally  .re  small  enough  for  their 
shape  to  be  strongly  influenced  by  surface  tension:  such  cavities 
usually  appear  in  the  fonn  of  spherical  or  nearly  spherical  bubbles. 
Cavitation  of  this  kind-sometimes  called  bulk  cavitation  because  of 
the  presence  of  a  large  number  of  bubbles  at  the  same  time-occurs 
when  the  pressure  in  a  large  volume  of  liquid  is  gradually  reduced  below 
the  liquid  vapor  pressure  or  the  saturation  pressure  of  dissolved  gases 
(effervescence);  or  when  a  tensile  stress  wave  propagate,  though  the 
liquid.  Depending  on  whether  the  ambient  pressure  field  leading  to 
cavitation  is  semistatic  or  dynamic,  the  cavitation  formation  and  growth 
process  is  governed  by  different  parameters.  A  slowly  varying  ambient 
pressure  field  usually  leads  to  the  appearance  of  gas,  or  gas-v.por- 
mixture  filled  bubbles,  while  large  magnitude  tensile  stress  waves 
produce  spalling  in  the  liquid,  which  manifests  the  formation  o,  nearly 
pure  vapor  bubbles.  Pure  vapor  cavitation,  unlike  gas  bubble  cavitation, 
cannot  be  eliminated  by  degassing  a  liquid. 

Pure  vapor  cavitation  is  of  particular  interest  in  „quid  explosive 
sensitivity  testing  because  large  magnitude  negative  stress  field,  are 
generally  associated  with  the  event  of  shock  impulses  propagating  through 
the  test  liquid.  The  negative  stress  fields  are  produced  indirectly  m 
the  test  liquid  as  a  consequence  of  primary  shock  propagation. 
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The  overall  shock  inpulse  consists  of  several  pulses.  First  to 
enter  the  test  liquid  is,  of  coursei  the  primary  shock.  This  primary 
peaked  compression  wave  is  subsequently  followed  by  secondary  and  higher 
order  waves  representing  the  first  and  higher  order  reflections  of  the 
various  attenuator  surfaces.  Cavitation  at  the  attenuator-liquid  inter¬ 
face  caused  by  rarefaction  waves  reflected  from  attenuator  free  surfaces 
has  been  observed  in  several  experiments.  2 > 9> 10 

As  the  main  shock  is  transmitted  to  the  liquid  under  testi  it 
produces  local  high  pressure  regions  in  the  liquid  which  in  turn  excite 
several  vibrational  modes  of  the  container  wall,  similar  to  the  familiar 
water-hammer  phenomena. 11  It  Is  assumed  here  that  the  stress  levels  set 
up  in  the  wall  do  not  exceed  the  elastic  limit  of  the  wall  material. 

The  wall  stress  waves  in  turn  cause  the  liquid-solid  boundary  to  oscillate 
and  radiate  pressure  waves  into  the  bulk  of  the  liquid,  giving  rise  to 
a  feed  back  loop  in  the  experimental  system. 

In  terms  of  magnitude  of  the  wall  deflections,  which  determines  the 
strength  of  the  wall-radiated  pressure  waves  in  the  liquid,  the  wall 
flexural  waves  and  wall  shear  waves  are  of  most  Importance.  If  the 
characteristic  wall  wave  speed  is  higher  than  the  speed  of  the  primary 
shock  in  the  liquid,  the  wall  vibrations  will  propagate  ahead  of  the 
primary  shock,  and  cause  precursor  pressure  wave  disturbances.  Positive 
and  negative  stress  levels  on  the  order  of  1/10  kbar  in  the  liquid  can 
be  produced  by  the  wall  precursor  wave  systems.  These  stress  levels 
are  sufficient  to  produce  cavitation  in  most  low  viscosity  liquids  such 
as  water. 

There  are  basically  4  different  sources  of  pure  vapor  cavitation 
in  liquid  explosive  sensitivity  testing.  In  addition  to  wall  precursor 
cavitation, 12  and  cavitation  existin';  at  the  attenuator  liquid  interface, 
cavitation  can  be  produced  in  weak  containers  due  to  radial,  plastic 
wall  deformation  behind  the  primary  shock13  and  also  in  the  bulk  of  the 
liquid,  as  the  primary,  peaked  shock  wave  is  reflected  off  the  free 
liquid  surface  cavitation  behind  the  primary  shock  should  start  at  th' 
wall  and  then  appear  along  the  Math  lines  toward  the  center  of  the 
container.  Precursor  cavitation  may  or  may  not  start  at  the  walls, 
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depending  on  the  magnitude  of  the  wall  deflection.,  and  ia  most  likely 
to  be  observable  near  the  center  of  the  container  where  wave  interaction 
gives  rise  to  maximum  .tress  levels.  Both  cavitation  at  the  attenuator 
liquid-interface  and  near  the  free  surface  of  the  liquid  should  appear 
first  through  the  central  region  of  the  container.  Thus  the  location 
of  the  first  appearance  of  cavitation  may  give  a  clue  as  to  the  type  of 
negative  stress  wave. 

Very  little  data  are  available  for  predicting  the  single  pulse 
cavitation  threshold  stress  levels.  Experiment,  in  the  field  of  acoustic 
cavitation14  indicate  the  cavitation  threshold  conditions  depend  on 
the  size  of  the  cavitation  nuclei  present  in  the  liquid,  on  the  pulse 
length,  15  on  the  liquid  viscosity,  and  other  factors. 14  Bull*  found 
from  single-pulse  experiments  in  various  liquids  that  the  negative 
cavitation  threshold  stresses  varied  with  viscosity  raised  to  the  1/5 
power.  However,  his  results  are  somewhat  in  question  since  he  did  not 
determine  the  size  of  the  cavitation  nuclei  present  in  his  test  liquids, 
and  it  is  the  latter  parameter  which  is  thought  to  have  a  major  effect 
on  the  cavitation  threshold. 

The  energy  necessary  for  initiation  of  the  exothermic  decomposition 
process  inside  the  vaporous  cavity  and/or  in  the  surrounding  liquid 
is  provided  by  the  bubble  collapse  process.  Bubble  collapse  is  brought 
about  by  either  overexpansion  of  the  bubble  beyond  the  diameter  in 
the  equilibrium  with  the  ambient  pressure,  or  by  a  second  compression 
pulse.  Cavities  collapsing  with  sufficient  violence  to  reach  spontaneous 
ignition  conditions  inside  the  bubble  can  constitute  hot  spots  in  the 
explosive  liquid.  In  the  vaporous  cavity,  the  effect  of  vapor  compression 
can  be  reduced  by  vapor  condensation  at  the  liquid-vapor  interface. 
However,  condensation  produces  a  large  amount  of  energy  release  which 
is  difficult  to  dissipate  in  an  environment  undergoing  shock  compression. 
The  amount  of  vapor  that  can  be  condensed  therefore  is  limited  by  the 
energy  losses  to  the  ambient  liquid.  In  general,  the  development  of 
hot  spots  in  a  liquid  explosive  can  have  three  possible  consequences: 
the  energy  losses  from  a  small,  iaolated  hot  spot  to  the  liquid  may  be 
too  large  to  sustain  a  reaction  and  the  hot  spot  may  disappear;  the 
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hot  spot  may  Initiate  slow  deflagration  burning  in  the  liquid,  which 
could  account  for  low-velocity  detonation  if  the  hot  spots  are  homo¬ 
geneously  distributed  across  the  container  c ross  section;  or,  under 
appropriate  experimental  conditions,  the  hot  spot  may  initiate  a  high- 
velocity  detonation. 

In  summary,  two  types  of  cavitation  are  of  interest  in  shock 
sensitivity  experiments:  gas  or  gas-vapor  filled  bubbles  and  pure 
vaporous  cavitation.  While  the  above  discussion  dealt  mostly  with 
pure  vaporous  cavitation,  the  presence  of  gas-filled  bubbles  could 
greatly  increase  the  shock  impulse  sensitivity  characteristic  of  the 
liquid  explosive.  However,  gas  bubbles  can  be  easily  eliminated  by 
degassing  the  liquid.  Vaporous  cavitation  visually  appears  in  the 
presence  of  peaked  compression  pulses  such  as  those  used  in  gap  sen¬ 
sitivity  testing.  The  occurrence  of  vaporous  cavitation  depends  not 
only  on  the  physical  parameters  of  the  liquid,  such  as  viscosity, 
surface  tension,  etc.,  but  also  on  the  degree  of  purity  of  the  liquids, 
the  size  of  heterogeneous  cavitation  sites  present,  and  the  negative 
stress  level  in  the  liquid.  Some  of  the  effects  of  the  physical  liquid 
parameters  on  the  sensitivity  will  be  discussed  in  later  sections. 

Vaporous  cavitation  can  give  rise  to  detonation  initiation  phenomena 
similar  to  the  gas  filled  bubbles  but,  in  addition,  the  observed  existence 
of  a  cavitation  front  velocity  is  intriguing  as  a  possible  support 
for  the  cavitation  hypothesis  of  LVD. 

C.  Experimental  Studies 

1.  Cavitation  and  Shock  Sensitivity  Testing 

As  stated  above,  we  have  been  investigating  shock  wave  interaction 
in  LVD  gap  sensitivity  testing.  Previous  reports  (2,3)  have  described 
some  of  these  efforts.  Recently  we  have  obtained  shadowgraphlc  pictures 
of  cylindrical  steel  containers  of  nitromethane  under  typical  gap  test 
conditions.  Plexiglas^  windows  2.54  x  0,63  cm  allowed  visual  pictures 
to  be  taken  of  the  liquid  under  shock  stress  conditions.  The  tubes 
were  2.54  cm  in  diameter  and  10.4  cm  in  length  with  wall  thickness  of 
0.63  cm.  The  results  of  these  studies  are  shown  in  Fig.  3,  where  eight 

^Trademark,  Rohm  &  Haas  Co. 
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FIGURE  3  CAVITATION  IN  SHOCKED  NITROMETHANE  (U) 
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selected  frames  from  a  twenty-five-frame  sequence  are  listed.  The 
initial  shock  was  produced  by  a  100-gram  tetryl  booster  and  attenuated 
by  2.54  cm  of  Plexiglas.  The  window  was  situated  midway  up  the  tube. 

The  framing  rate  was  2.08  ^sec/frame  and  the  number  on  each  frame 
indicates  the  number  of  nsec  from  initiation.  The  shock  wave  arrives 
at  the  Plexiglas-tube  interface  at  19  nsec.  Crossed  precursor  waves 
are  seen  in  the  30.4  nsec  frame.  (The  photographs  reproduced  here 
are  half-tone  prints  which  represent  a  reduction  from  a  1/128  gray 
scale  for  the  negative  to  1/12  for  these  prints.  Therefore  some  of  the 
details  seen  in  the  negatives  are  lost  in  the  reproduction.)  Two 
frames  later  (34.6  nsec)  bubbles  are  clearly  visible  following  the  passage 
of  the  precursor  waves  (seen  at  top  of  window).  The  primary  liquid 
shock  arrives  in  the  fifth  frame  (42.9  nsec)  and  the  bubbles  are  com¬ 
pressed  by  its  passage;  howeveri  not  all  of  the  bubbles  appear  to  be 
compressed  by  the  liquid  shock.  In  the  later  frames  the  Plexiglas 
windows  show  signs  of  cracking. 

The  precursor  wave  has  an  average  velocity  of  about  2.5  mm/nsec; 
however  one  would  expect  this  velocity  to  be  nearer  that  of  the  steel 
wall  sonic  velocity  of  5.0  mm/nsec.  The  bubble  front  velocity  varies 
from  2.0  to  3.5  mm/nsec  and  appears  to  be  Increasing  with  time.  The 
liquid  shock  velocity,  as  expected,  decreases  from  2.0  to  ~  1.0  mm/nsec 
during  its  passage  across  the  viewing  port.  These  preliminary  results 
are  in  qualitative  agreement  with  that  outlined  above  concerning 
cavitation  in  gap  test  experiments.  More  quantitative  work  is  needed, 
however,  before  a  detailed  picture  of  cavitation  can  be  obtained  and 
compared  to  theoretical  studies.  We  are  planning  to  refine  our  techniques 
to  include  Schllern  photography  which  is  more  sensitive  to  the  density 
gradient  in  a  stress  wave. 

2.  Effects  of  Confinement  on  LVD  Gap  Sensitivity 

In  the  last  report9  we  presented  preliminary  data  on  the  effects 
of  wall  thickness  on  the  LVD  gap  sensitivity.  The  analyzed  results 
(Table  2)  show  that  the  apparent  LVD  gap  sensitivity  of  1,2-DP  increases 
with  decreasing  wall  thickness.  However,  there  does  appear  to  be  a 
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Table  2 


LVD  GAP  SENSITIVITY  OF  1,2- DP:  WALL  THICKNESS  STUDIES 
(STEEL  TUBES  1.27  cm  i.d. ,  10-ea  LENGTH) 


Wall 

LVD  Gap 

LVD  Transit  Velocity 
- (mm/asec) 

(cm) 

(cm  Plexiglas) 

0.63 

30  <  -  >  40.6 

2.3, 

0.23 

61  <  -  >  122 

2.4, 

0. 12 

91.5  <  -  >  107 

LVD  unstable 

m  2.4 

minimum  thlckneee  below  which  the  LVD  becoaee  unstable  and  doea  not 
propagate.  These  facts  are  further  evidence  of  the  coupling  of  the 
LVD  reaction  wave  with  the  confinement  walla.  It  appears  that  if  the 
walls  become  too  thin  they  cannot  support  a  stable  LVD  wave.  Lateral 
wall  movement  nay  occur  so  rapidly  that  the  reactions  are  quenched 
before  the  energy  Is  fully  released.  Figure  4  depicts  four  selected 
frames  from  a  high  speed  framing  camera  sequence  of  an  LVD  gap  sensitivity 
test  using  1,2-DP.  The  tuba  wall  thickness  is  0.12  cm.  Unstable  LVD 
is  seen  to  begin  in  the  535  psec  frame  as  the  tube  begins  to  swell. 

By  552  psec  the  reaction  has  proceeded  unstably  to  the  top  of  the  tube 
and  also  has  split  the  tube  walls.  In  the  thicker-walled  tests  the 
reaction  always  began  at  the  bottom  and  proceeded  normally  to  the  top 
of  the  tube.  The  initiation,  however,  was  delayed  on  the  order  of 
100  y,sec  after  the  initial  donor  shock  wave  had  reached  the  tube- 
attenuator  Interface.  This  delay  seems  to  be  characteristic  of  the 
LVD  reaction  in  1,2-DP  when  large  gaps  (>  30  cm)  are  used. 

We  have  completed  some  of  the  preliminary  studies  on  the  effect 
of  charge  diameter  on  the  LVD  gap  sensitivity.  Using  IBA  as  a  model 
compound,  we  determined  the  LVD  gap  sensitivity  in  2.54  and  1.27  cm 
diameter  steel  tubes  (wall  thickness  0.22  cm).  The  LVD  gap  sensitivity 
appears  to  Increase  with  increasing  charge  diameter.  However  these 
results  may  only  reflect  the  change  in  strength  of  the  wall  as  the 
diameter  is  increased.  We  are  planning  to  test  this  idea  by  a  set  of 
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FIGURE  4  UNSTABLE  LVD  IN  1,2-DP  (0.12  cm  tubes)  (U) 
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experiments  in  which  the  ratio  of  the  inside  diameter  (D)  to  wall  thickness 
(t)  remains  constant.  We  will  do  this  for  two  values  of  D/t.  These 
results  will  be  discussed  in  the  next  report. 

3-  Study  of  Tracer  Dye  Technique 

In  conjunction  with  our  studies  of  the  shock  wave  interaction 
occurring  in  gap  sensitivity  testing,  we  investigated  the  possibility 
of  using  a  tracer  dye  technique  to  follow  the  flow  field  behind  these 
shock  interactions.  This  technique  was  developed  by  Popovich  and 
Hummel 17  to  follow  developing  boundary  layers.  The  tracer  material, 

2-(2' ,4 '-dinitrobenzyl) -pyridine  (DNBP) ,  is  colorless  until  it  is  activated 
to  a  colored  tautomeric  form  by  uv  light.  Since  the  colored  form  will 
persist  for  several  milliseconds,  this  material  can  be  used  as  a  tracer 
in  a  flow  field.  We  had  planned  to  flash-photolyze  a  thin-colored 
line  before  the  arrival  of  the  weak  shock  waves  and  then  follow  the 
shock  flow  field  by  high  speed  photography.  The  colorless  material 
strongly  absorbs  only  about  3000  k,  but  nitro  and  nitrate  groups 
common  to  most  liquid  explosives  strongly  absorb  in  this  region  (See 
Fig.  5).  Therefore,  this  technique  would  require  the  use  of  a  non¬ 
explosive  liquid.  The  difluoramlno  compound  could  be  used  since  they 
do  not  begin  to  strongly  absorb  until  2000  A,  but  they  may  not  be  com¬ 
patible  with  the  dye  material.  In  addition,  as  Fig.  6  shows,  window 
material  would  be  limited  to  uv  transparent  material  like  Vycor.t 
Plexiglas  absorbs  too  strongly  in  the  3000  k  region  to  be  useful.  The 
colored  fonn  of  DNBP  (li„e  4,  Fig.  6)  has  its  absorption  centered  about 
6000  k  so  it  could  be  photographed. 

Using  glass  material  for  windows  does  present  some  difficult  problems. 
Hie  glass  fracture  rate  under  this  condition  is  very  high  (5  mm/^sec) 
and  this  fracturing  would  likely  interfere  with  subsequent  pictures. 

The  window  material  requirement  presents  a  form!  able  barrier  to  the 
development  of  DNBP  as  a  tracer  dye  for  shock  experiments  with  explosive 
liquids.  On  this  basis  we  decided  to  abandon  any  further  development 
of  this  technique. 


^Corning  Glass  Works. 
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FIGURE  5  A8SORPTION  SPECTRA  OF  THREE  COMPOUNDS 
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FIGURE  6  ABSORPTION  SPECTRA  OF  WINDOW  MATERIAL  AND  DNBP 
(Colored  and  Colorless) 
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D. 


Discussion 


Our  studies  of  LVD  sensitivity  in  the  dif luoramino  compounds 
have  been  guided  by  an  attempt  to  discover  the  determining  factors. 

Using  the  gap  sensitivity  test  we  have  measured  the  effect  of  confinement 
on  LVD  sensitivity.  These  effects  are  summarized  In  Table  3. 

Table  3 

EFFECT  OF  CONFINEMENT  ON  THE  LVD  GAP  SENSITIVITY 
OF  1,2-DP 

Confinement  Parameter  _ Effect _ 

Confinement  geometry1*  Increased  sensitivity 

in  cylindrical  geometry. 

Sonic  Velocity1*  Increased  sensitivity 

with  increased  sonic 
velocity. 

Wall  thickness*  Increased  sensitivity 

with  decreased  wall 
thickness  (unstable  in 
very  thin  walls) 

It  should  be  understood  that  these  confinement  effects  may  not 
generalize  to  all  LVD-sustalnlng  liquids,  but  we  hope  they  can  act  as 
a  useful  guide  to  a  probable  mechanism.  They  support  the  proposition 
that  in  LVD  there  is  a  coupling  of  the  explosive  and  its  confinement. 

If  this  coupling  is  accomplished  by  precursor  waves,  then  increasing 
the  ratio  of  the  wall  sonic  velocity  to  that  of  the  liquid  would 
magnify  the  effects  of  these  waves.  Decreasing  the  wall  thickness 
(strength)  also  would  act  to  magnify  these  effects.  This  is  of  course 
what  is  found  experimentally.  In  cylindrical  geometry  there  would  be 
a  stronger  focusing  effect  on  the  precursor  waves  than  in  rectangular 
confinement.  We  have  noted  an  increase  in  the  LVD  gap  sensitivity  in 
cylindrical  geometry.  The  charge  diameter  studies  are  not  yet  complete 
so  we  will  refrain  from  comment  on  their  meaning  at  this  time. 
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We  have  also  studied  the  effects  of  chemical  properties  on  LVD 
gap  sensitivity.  By  comparing  results  obtained  for  various  difluoramino 
isomers,  we  were  able  to  show  that  compounds  with  the  fastest  low 
pressure  reaction  rate  were  the  most  sensitive  to  LVD.*  In  addition, 
we  compared  these  results  with  other  sensitivity  tests.  These  comparisons 
(Table  4)  show  that  there  appears  to  exist  a  correlation  between  the 
gas  phase  reaction  rate  and  the  LVD,  drop  weight,  and  spark  sensitivity. 

One  conclusion  that  can  be  drawn  from  this  correlation  is  that  the 
reactions  that  lead  to  positive  results  in  these  tests  are  similar. 

The  reaction  of  importance  under  these  conditions  is  the  scission  of  the 
C-N  bond.  This  reaction  is  favored  at  low  pressure  end  lower  temperatures. 


Table  4 


CORRELATION  OF  SENSITIVITY  TESTS  WITH  REACTION  RATES 


Gas-Phase 
[  Reaction  Rate 

Sensitivity  Test 

Reaction  Time 
HVD  Gap 

LVD  Gap 

Drop  Weight 

Spark 

2,2-DP  , 

a 

2,2-DP 

2,2-DP 

2,2-DP 

2,2-DP 

b 

1,2-DP 

1,2-DP 

1 , 2-DP 

1,2-DP 

1,2-DP 

1,3-DP 

1,3-DP 

— 

1 

t  =  increasing  sensitivity  or  reactivity 
1  =  increasing  reactivity. 


On  the  other  hand  the  high  velocity  detonation  (HVD)  reaction  time  corre¬ 
lates  very  well  with  the  rate  of  elimination  of  hydrogen  fluoride  (HF). 
This  reaction  is  favored  under  the  conditions  of  high  pressure  and 
temperature  that  prevail  in  HVD. 3 

If  the  initiation  of  the  LVD  reaction  is  a  low  pressure-lower 
temperature  phenomenon,  what  are  the  conditions  that  can  provide  this 
environment?  In  Section  B-2  we  discussed  some  of  the  parameters  that 

fThe  chemical  kinetic  data  were  obtained  by  D.  S.  Ross  as  a  part  of  our 
overt.ll  studies  of  sensitivity  fundamentals.  See  refs.  2  and  3  and 
Section  IV  for  the  details  of  these  studies. 
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lead  to  cavitation.  Precursor  waves  can  lead  to  cavitation  ahead  of 
the  primary  liquid  shock.  The  interaction  between  the  liquid  shock 
and  cavitation  could  lead  to  initiation  of  LVD.  However,  we  have  noted 
in  many  LVD  tests  that  there  are  long  delays  between  the  initial  shock 
arrival  and  initiation  of  LVD.  Therefore  the  situation  is  more  complicated 
than  this  simple  explanation  would  indicate. 


If,  however,  cavitation  is  an  important  part  of  LVD  then  certain 
physical  characteristics  of  the  explosive  should  be  the  determining 
factors  of  LVD  initiation.  These  factors  include  the  heat  release 
(AH  ),  the  heat  vaporization  (AH^),  the  surface  tension,  the  heat 
capacity,  the  sonic  velocity,  the  viscosity,  the  gas  solubility,  and 
the  chemical  reaction  kinetics.  It  appears  that  two  of  the  more  important 
parameters  in  LVD  would  be  the  heat  of  reaction  and  the  heat  of  vaporization. 
AH  would  supply  the  energy  to  drive  the  LVD  wave  and  AHy  would  remove 
energy  in  the  cavitation  process.  Table  5  lists  several  compounds  for 
which  we  have  calculated  the  heat  of  reaction  based  on  the  assumption 
of  simple  reaction  products.  The  values  of  AHy  in  Table  5  were  obtained 
from  the  literature  or  were  calculated  from  the  vapor  pressure  curves. 

We  have  also  included  the  LVD  gap  sensitivities  of  these  selected 
compounds . 


Figure  7  is  a  plot  of  AHr/AHv  as  a  function  of  the  LVD  gap  sensitivity 

for  the  compounds  in  Table  5.  We  have  included  the  two  dinit ropropane 

isomers  as  we  are  planning  to  obtain  their  LVD  gap  sensitivity.  The 

correlation  between  AH  /AH  and  the  LVD  gap  sensitivity  is  very  good 
r  v 

considering  this  simple  approach.  This  would  further  indicate  that 
perhaps  these  physical  chemical  factors  are  the  dominant  ones  in  LVD. 
However,  as  we  have  stated  above  the  reaction  rate  also  plays  an 
important  role  in  LVD  initiation.  We  have  included  nitroglycerine  in 
Table  5  but  not  in  Fig.  7  because  we  do  not  have  a  value  for  its  LVD 
gap  sensitivity.  But  we  do  know  that  its  LVD  sensitivity  value  is 
very  large  which  would  be  indicated  by  the  rather  high  value  for  AHr/AHv- 

The  purpose  of  correlation  is  to  predict  values  for  new  materials. 

So  we  have  estimated  the  LVD  gap  sensitivity  of  the  dinit ropropanes 
by  calculating  their  AHr/AHv>  Based  on  this  estimation  we  would  conclude 
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Table  5 


SOME  THERMODYNAMIC  PROPERTIES  AND  TOE  LVD  GAP 
SENSITIVITIES  OF  SELECTED  EXPLOSIVES 


Compound 

Heat  of 
Formation 

AHf 

Heat  of 
Vaporization 
AHV 

Heat  of 
Reaction 
AHr 

LVD  Gap 

Sensit ivitya 

Kcal/mole 

Kcal /mole 

Kcal /mole 

cm  Plexiglas 

Nltromethane 

-27.03 

9.09 

-69.38 

0 

Ethyl  Nitrate 

-44.3 

8.  7 

-100.2 

3. 8-5.1 

TOETN 

-105.8 

20.7 

-273.2 

5-10 

1 , 2-DP 

-48.0 

I  I 

-205 

30-40 

1 , 1-Dinit  ropropane 

-40. 78 

-152.04 

b 

1-5  Est . 

2, 2-Dinit ropropane 

-44.87 

-147.95 

5-10  Est . b 

Nitroglycerine 

-90.75 

-1214.5 

very  large 

a 

Data  obtained  using  either  camera  or  probes  as  detectors, 
b 

Values  estimated  from  Fig.  7. 


that  these  are  not  too  sensitive.  In  fact>  the  1,1-isomer  may  not  show 
any  LVD  at  all.  ’Ye  are  planning  to  test  them  for  LVD  as  soon  as  possible 
to  see  if  we  have  made  a  correct  estimate.  It  should  be  emphasized 
again  that  the  parameter  AH^/AH^  is  not  the  only  one  that  may  be  important. 
In  fact  the  chemical  reaction  kinetics  should  be  taken  into  account. 
Nltromethane,  which  does  not  undergo  LVD,  has  several  of  these  factors 
against  it.  Its  heat  of  reaction  is  low.  Although  its  AHv  is  reasonable, 
its  reaction  kinetics  are  known  to  be  very  slow. 

We  hope  in  the  future  to  develop  a  parameter  that  include  many 
other  factors  of  importance  in  LVD  initiation.  In  this  way  one  may  be 
able  to  calculate,  a  priori,  the  LVD  gap  sensitivity  of  a  new  material 
from  knowledge  of  its  physical  and  chemical  properties. 
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E, 


Future  Work 


We  will  continue  our  efforts  to  determine  the  important  parameters 
in  LVD  gap  sensitivity.  From  these  studies  we  hope  to  derive  more 
information  about  low-velocity  detonation  in  general  as  well  as  LVD 
gap  sensitivity  in  particular.  During  the  next  period  we  will  continue 
theoretical  studies  of  cavitation  effects  as  well  as  experimental  programs 
to  complement  these  studies. 

We  will  continue  our  efforts  to  obtain  charge  diameter  effects  on 
LVD  gap  sensitivity.  In  addition,  we  will  study  the  effects  of  charge 
length  on  LVD  sensitivity.  The  materials  are  now  on  hand  for  these 
st udies. 

Some  effort  will  be  expended  in  LVD  studies  of  the  dlnitropropanes. 

We  are  particularly  Intrigued  to  discover  if  our  prediction  of  their 
LVD  gap  sensitivities  are  correct. 

We  will  also  attempt  to  further  develop  correlation  parameters  for 
LVD  gap  sensitivity.  In  particular,  we  would  like  to  Include  reaction 
kinetics  and  attempt  to  discovery  how  other  parameters  correlate  with 
LVD  gap  sensitivity;  the  ultimate  aim  of  such  studies  being  the  pre¬ 
diction  a  priori  of  the  sensitivity  of  new  materials. 
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Ill  PHYSICS  AND  CHEMISTRY  OF  DETONATION 


(M.  Cowperthwaite  and  R.  Shaw) 


A.  Introduction 

This  program  is  concerned  with  the  transient  detonation  phenomena, 
shock  initiation  and  failure.  During  the  present  reporting  period  we 
have  concentrated  on  shock  initiation.  We  have  devoted  most  effort 
to  solving  one  of  the  major  problems  in  shock  initiation-calculation 
of  the  shock  temperature.  Calculation  of  shock  temperature  is  important 
in  shock-init iation  studies  of  liquid  explosives  because  pressure- 
volume-temperature  (p-v-T)  equations  of  state  of  liquids  in  the  kilobar 
region  are  not  known.  Use  of  the  Walsh-Christian  method1  developed 
for  metals  is  obviously  limited  by  the  assumptions  upon  which  it  is 
based.  The  assumption  of  constant  specific  heat  at  constant  volume  C 

v 

is  not  valid  for  molecular  liquids  with  internal  degrees  of  freedom, 

neither  is  the  assumption  that  numerical  differences  between  C  and 

v 

specific  heat  at  constant  pressure  Cp  can  be  neglected  for  calculational 

purposes.  It  is  necessary  to  account  for  those  differences  between 

liquids  and  solids  to  make  the  Walsh-Christian  method  more  applicable 

to  liquids.  The  present  work  provides  a  new  method  of  calculating 

shock  temperatures.  It  is  based  on  a  variational  analysis  which  shows 

that  calculated  shock  temperature  is  much  more  sensitive  to  changes 

in  Cv  than  to  changes  inOp/ST)^.  Shock  temperatures  for  various 

liquids  were  calculated  with  both  constant  and  variable  C  and  those 

v 

for  carbon  tetrachloride  were  compared  with  brightness  temperatures 
measured  by  other  workers.2)3 

B.  The  Walsh-Christian  Method  of  Calculating  Shock  Temperature 

Cowperthwaite4  has  discussed  the  significance  of  two  equations 
of  state,  the  Walsh-Christian  method  (constant  Cy)  used  for  metals1 
and  the  Rice-Walsh  method  (constant  C  )  used  for  water,5  For  the 
present  work  the  Walsh-Christian  method  was  preferred  because  it 
required  less  high  pressure  data. 
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Briefly,  the  basis  of  the  method  is  as  follows.  From  the  Rankine- 
Hugoniot  relations, 

de  =  1/2  (p  -  p)dv  +  1/2  (v  -  v)  dp  (1) 

'  o  o 

where  e  is  the  internal  energy,  pQ  is  the  initial  pressure  (assumed  to 

be  negligible  compared  to  the  pressure  at  p),  and  v  is  the  initial 

o 

volume. 


Further, 


de  =  1/2  pdv  +  1/2  (v  -  v)  dp 


de  =  (3e/3T)vdT  +  (de/^v)  dv 


(21 


(3) 


where  T  is  the  temperature.  Substituting  C^,  the  heat  capacity  at 
constant  volume  for  Oe/3T)v,  and  [T(3p/5T>v  -  p]  for  Oe/^v)^  in 
Eq.  (3), 

de  =  CydT  +  f T(3p/$T) ^  -  p]dv  (4) 

Combining  Eqs.  (2)  and  (4), 

CydT  +  T(3p/3T)vdv  =  l/2pdv  +  1/2(vq  -  v)dp  (5) 

dividing  throughout  by  Cvdv,  we  obtain  the  differential  equation  along 
the  Hugoniot 


dT  T($p/ST)V  =  _1_  _  v)  dfi] 

dv  C  2C  ’’  o  dvJ 

v  v 


(6) 


Walsh  and  Christian  showed  that  ±t  (^p/dT>v  and  C^  are  constant 
Eq.  (6)  can  be  solved  explicitly  to  give 


b(v  -  v)  -bv 

T  =  T  e  0  +  J  ebVF(v)dv 

o  2C  v 


(7) 


where  b  =  Op/&T>v/C  and  F(v)  =  p  +  (vq  -  v)(dp/dv). 
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A  BASIC  computer  program  was  written  by  B.  Y.  Lew  of  SRI  to  calculate 
shock  temperatures  using  Eq.  (7).  In  the  program  the  integral  was 
evaluated  by  a  trapezoidal  approximation.  The  program  was  checked  as 
follows.  Walsh  and  Christian's  input  data  for  b,  C  and  F(v)  for 
copper  were  used  and  excellent  agreement  was  obtained  with  their  shock 
temperatures.  In  addition,  C.  L.  Mader  of  Los  Alamos  Scientific 
Laboratory  kindly  used  our  input  data  for  nitromethane  with  his  program 
and  obtained  agreement  with  our  results. 

Another  way  of  calculating  the  shock  temperature  is  to  integrate 
the  differential  Eq.  (6)  using  a  Runge-Kutta  technique.  This  is  a  more 
generally  useful  method  because  any  function  of  C  and  (dp/3T)  can 
be  used.  Of  course,  if  constant  values  of  Cy  and  Op/3T>v  are  used, 
then  the  shock  temperatures  calculated  by  the  Runge-Kutta  integration 
of  Eq.  (6)  will  be  the  same  as  those  calculated  by  the  trapezoidal 
evaluation  of  the  integral  in  Eq.  (7).  This  agreement  was  confirmed 
by  the  calculations  for  carbon  tetrachloride  shown  in  Table  6.  The 
input  data  is  in  Table  7.  The  p-v  Hugoniot  from  which  F(v)  is  derived 
is  obtained  by  using  the  Rankine-Hugoniot  relations  to  transform  the 
universal  liquid  liu„oniut 


Us  =  u>Co  +  u*up  (8) 

where  Ug  is  the  shock  velocity,  Up  is  the  particle  velocity,  cq  is  the 
sound  velocity  at  ambient,  and  Uj  and  u2  are  constants. 

In  Fig.  8  the  shock  temperature  for  carbon  tetrachloride  calculated 
using  constant  Cy  and  Op/aT^  is  compared  with  the  "brightness"  tempera¬ 
tures,  previously  measured  by  Voskoboinikov  and  Bogomolov2  and  by 
Ramsay. 3 

Clearly,  the  experimental  and  calculated  results  are  in  poor 
agreement.  We  consider  it  likely  that  the  calculated  values  are 

"too  high"  or,  in  other  words,  the  model  for  the  calculated  temperatures 
is  too  simple. 
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TEMPERATURE  —  °K 


FIGURE  8  SHOCK  TEMPERATURE  FOR  CARBON  TETRACHLORIDE.  Comparison  of  the 
calculated  values  with  those  obtained  experimentally  bhthe  "brightness"  method. 
The  line  Cv  was  calculated  in  the  present  work  using  constant  Cv.  The  line  MCp 
was  calculated  by  Mader  (Ref.  3)  using  constant  Cy  but  using  Cp  for  the  value 
of  Cy.  The  dashed  line  CV(T)  was  calculated  in  the  present  work  using  Cy  as  a 
function  of  temperature.  The  input  data  for  the  calculation  are  in  Table  2. 


32 


C-  Variational  Analysis  of  the  Walsh-Chrlstlan  Method  (Sensitivity 

°^8  Calculated  Temperatures  to  tha  Values  CST»F  for  ^p7§f^:, 

Of  the'  three  variable.  (ap/aT)y,  Cy,  and  F(v)  In  the  differential 
Eq.  (6),  F(v)  la  probably  the  beat  eatabllahed.  It  May  be  assumed  to 
be  known.  (It  will  later  be  shown  empirically  that  the  calculated 
temperature  as  a  function  of  shock  pressure  Is  almost  unaffected  by 
using  different  fonaa  for  F(v),) 

Concerning  (bp/bT)v  and  Cv,  consider  integration  between  v0  and 
V!.  At  v,,  let  Tg  be  the  temperature  on  the  Hugoniot  and  let  T  be  the 
temperature  on  the  isentrope,  given  by  T,  =  IQ  exp  b(vQ  -  Vl).  'consider 
the  shock  temperature  at  v,  to  be  a  function  of  b  and  C  ,  then  f-.-om 
Eq.  (7)  v 

-bv,  Vl 

Th(v,,  b,  Cy)  =  Tg  +  — -  J*  ebVF(v)dv  (9) 


-  t.  *35-  ;  •b<v-v‘> 


F(v)dv  (10) 


The  slope  of  the  Tg  versus  (ap/3T)  curve  Is  given 


3  0p/^t)v  ~  Ts  <vo  Vl>  ^ Op/aT)..  +  ac.  /  (v_vi>  C3(3p^T)  ]  ®b(v  Vl)F(v) 


v  v 


now  b  =  (ap/3T)v/cy  .  .  ab/a(ap/aT)y  =  i/c 

and  5TH  VVV>>  I 

and  S(ap/av)y  =  — —  +  ^7 


where  I  =  J  (v-vt)  eb(v-v‘) 


F(v)dv 
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The  RHS  of  Eq.  (12)  is  positive,  so  increasing  (3p/3T>v  will 

increase  T„.  Similarly,  the  slope  of  the  T  versus  C  wave  is  given 
H  H  V 

by 


ac 


Wv‘>  H 


1 

2C 


1  (*  ,  ,  Sb  Mv-v,)-,.  .  . 

“  J  (v'Vl>  3c.  6  F(v)dv 


P*  b(v-vt) 

e  1  F'v)dv 


(13) 


Now  b  =  Op^T)v/Cv 

ayacv  = 


ab/ac  =  -  Op/st)  /c  2  =  -  b/c 


-  b 


T  (v  -Vj) 
so1 

i 

1 

*-3 

in 

C 

+  2C  2 

c  l 

v 

V 

L  v  J 

9T„ 


( T  -  T  ) 
H  s 


a(ap/5T) 


(14) 


The  RHS  of  Eq,  (14)  is  negative  so  that  decreasing  results  in  an 
increase  in  T, 


Let  R=  OTH/3Cv)  /  (aTH/a(dp/ST)v) 


Then  R  = 


f  3Th  \ 

ft  - 

a(ap/at)v/ 

l  CV  / 

v/ 

^TiF 


a(ap/aT)v 


TH  '  TS  \ 

cvoya(aT)v  j 

T„  -  T  \ 

_ H  s _ j 

Ts(vo”Vl>  +  v(I/2Cv)  / 


(15) 
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For  a  lot  increase  in  Op/aT)^  i.e.  AOp/ST^  =  .l(aP/3T>v  and  a  10t 
decrease  in  Cy  i.e.  AC^  =  -.1  Cy  neglecting  second  and  higher 
derivatives. 


AT(C  ) 

_ _v _ _ 

ATOp/3T)v 


1 


(16) 


As  Th  is  always  greater  than  Ty  the  shock  temperature  is  more  sensitive 
to  Cy  than  to  0)p/aT)v.  Equations  (15)  and  (16)  have  been  evaluated 
for  carbon  tetrachloride  up  to  250  kbar.  The  results  in  Table  8  show 
that,  as  the  pressure  increases  up  the  Hugoniot,  C  becomes  increasingly 
more  important  than  Op/&T)  . 

The  above  conclusions  were  confirmed  empirically  as  shown  by  Fig.  9 
and  the  last  column  of  Table  8. 

To  test  the  effect  of  changing  the  p-v  Hugoniot,  the  shock  temperature 
was  calculated  for  four  different  "universal"  Hugoniots  for  carbon 
tet rachloride,  namely 

Us  =  U2Co  +  ^’"p  (ref-  2)'  us  =  l-25co  +  1. 7Up  (this  work,  arbitrary 
variation  from  1.2  to  1.25),  Ug  =  1.199c  +  1.672u  (ref.  7),  and 
U  =  1.31c  +  1. 61u  (ref.  8).  P 

s  O  P 

The  last  Hugoniot  had  previously  been  suggested8  as  a  better  "universal” 
fit  for  liquids.  The  results  in  Fig.  10  show  that  the  temperature 
as  a  function  of  shock  pressure  is  little  changed,  although  Table  9 
shows  that  the  p-v  Hugoniots  are  significantly  different. 

D.  Effect  of  Increasing  Pressure  and  Temperature  on  Op/3T)  and  C 

For  carbon  tetrachloride  at  atmospheric  pressure  around  room 
temperature  (3p/3T)v  increases  with  temperature.  9  But  from  Bridgman's 
work  on  the  p-v-T  relation  of  liquids  at  high  pressures,  Op/3T) 
decreases  with  increasing  temperature  at  constant  volume.  Bridgman 
made  a  detailed  investigation  of  the  previously  advanced  hypothesis 
that  (3p/3T)v  is  a  function  of  volume  only.  Although  the  relation 
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TEMPERATURE  —  °K 


PRESSURE  —  kbar 

TA-4501-4J2 

FIGURE  9  SHOCK  TEMPERATURE  OF  CARBON  TETRACHLORIDE  CALCULATED 
USING  THE  CONSTANT  Cv  MODEL.  Sensitivity  of  the  calculated 
temperatures  to  the  values  used  for  Cv  and  (dp/8T)v.  The  input  data 
are  in  Table  2. 
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CARBON  TETRACHLORIDE 

TEMPERATURE  —  °K 


PRESSURE  —  kbar 

T  A-410 1-431 


FIGURE  10  SHOCK  TEMPERATURE  FOR  CARBON  TETRACHLORIDE.  Sensitivity  of  the 
calculated  temperatures  to  the  form  of  the  Hugoniot.  The  input  data  are  In 
Table  2. 
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starts  to  break  down  at  high  pressures,  it  holds  well  at  low  and 
moderate  pressures  and  is  a  reasonable  approximation. 

If  Op/3T>v  is  a  function  of  volume  only,  then  it  follows  that 
is  a  function  of  temperature  only.  At  atmospheric  pressure  and 
around  room  temperature  for  carbon  tetrachloride  and  similar  liquids, 

Cv  increases  with  temperature.  Bridgman10  showed  that  at  temperatures 
up  to  about  450°K,  of  18  liquids  did  not  vary  much  with  pressure 
up  to  100  kbar. 

E.  The  CV(T)  Model 

From  the  preceding  sections,  two  main  conclusions  emerge  concerning 

the  shock  temperature  calculated  by  the  Walsh-Christian  method.  First, 

the  calculated  temperature  is  more  sensitive  to  variation  of  C  than  of 

v 

Op/dT)^.  Second,  as  the  temperature  and  pressure  increase  along  the 
Hugoniot,  it  seems  likely  that  (op/3T>v  will  at  first  increase,  then 
will  decrease,  whereas  it  is  probable  that  will  increase.  In  attempting 
to  develop  a  better  equation  of  state,  we  have  chosen  to  retain  the 
Walsh-Christian  assumption  that  (5p/5T>v  is  constant,  because  the  shock 
temperature  is  less  sensitive  to  (jp/dT>v  than  to  C^,  because  the  increase 
and  subsequent  decrease  of  Op/3T>v  will  tend  to  cancel  and,  most 
important  of  all,  we  have  little  intuition  as  to  the  variation  of  ( 5 p/3 T )  . 

On  the  other  hand,  the  calculated  shock  temperature  is  more  sensitive 
to  C^  than  (Sp/dT)^  and,  as  it  is  likely  that  will  increase  as  the 
temperature  and  pressure  are  increased,  it  seems  worthwhile  to  make  some 
assumptions  concerning  the  variation  of  with  temperature  and  pressure. 
From  the  discussion  in  the  previous  section  it  is  a  reasonable  assumption 
that  will  be  a  function  of  temperature  only.  Consider  the  heat 
capacity  of  an  ideal  gas.  For  polyatomic  molecules,  there  is  usually 
a  big  increase  in  as  the  temperature  is  raised  from  room  temperature 
to  several  thousand  degrees  Kelvin.  The  increase  in  Cy  is  due  primarily 
to  an  increase  in  the  vibrational  heat  capacity  as  the  oscillators 
become  more  classical  as  l-(hv/kT)  approaches  unity.  It  has  previously 
been  assumed  that  the  vibrational  heat  capacity  has  the  same  value  in 
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the  liquid  as  in  the  ideal 


gas. 


Therefore  we  have  assumed  that 


1 1 


CV(T)4  =  cv(298)t  + 


where  Cv(T)^  is  C 

298°K  and  AC0  i 
vg 

Because  C  and  C 

„  v  p 
AC0  =  AC0  . 
vg  pg 


v  of  the  Uduid  at  T.  Cv(298)^  is  Cy  of  the  liquid  at 
s  the  increase  in  C y  of  the  ideal  gas  from  298°K  to  T. 
of  an  ideal  gas  differ  only  by  R,  it  follows  that 


A  computer  program  was  written  to  integrate  the  differential 
Eq.  (7)  by  a  Runge-Kutta  technique.  If  Cy  and  (ap/&T)y  are  constant, 
then  the  Runge-Kutta  method  should  give  the  same  result"  as  the  trapezoidal 
integration  of  Eq.  (8).  This  was  checked  for  carbon  tetrachloride, 
the  temperatures  calculated  by  the  two  methods  agreeing  to  better  than 
1". 


F-  Results  of  Temperature  Calculations 

Shock  temperatures  were  calculated  for  carbon  tetrachloride, 
ni t romethane,  and  water. 


!*  Carbon  Tetrachloride 

The  shock  temperature  of  carbon  tetrachloride  was  calculated  using 
Cv(T>  and  other  input  data  as  described  in  Table  7.  The  results,  shown 
in  Fig.  8,  show  much  better  agreement  with  the  experimental  measurements 
than  do  the  temperatures  calculated  using  the  constant  value  of  C  . 

Mader12  obtained  good  agreement  with  the  experimental  results  usi^g  the 
Walsh-Christian  method  (constant  Cy) .  However,  the  agreement  is  fortuitous 
because  he  used  the  value  of  Cp  for  Cy  (0.198  cal  mole'2  deg-  lnstead  of 
0.141  cal  mole-1  deg-1). 


It  is  interesting  to  note  that  the  experimentally  observed  temperatures 
start  to  diverge  from  those  calculated  using  Cm  at  pressures  above 
about  150  kbar.  This  is  the  region  in  which  Dick6  observed  a  break  in 
the  p-v  Hugoniot  and  is  also  the  region  where  Mader13  calculated  that 
significant  amounts  of  decomposition  of  CC1.  into  C2C16  and  Cl2  take 
place. 
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2. 


Nitromethane 


Figure  11  shows  the  temperatures  calculated  for  nitromethane  using 

constant  C  and  C  (T).  The  input  data  are  in  Table  7.  The  results 
v  v 

are  compared  with  those  calculated  by  Enig  and  Petrone14  using  another 

equation  of  state.  From  their  equation  of  state,  Enig  and  Petrone 

obtained  C  as  a  function  of  temperature.  The  value  of  C  at  298°K 
v  v 

that  they  obtained  (0.24  cal  g~ 1  deg'1)  is  significantly  different  from 
the  literature-  value  (0.29  cal  g-1  deg-1).  More  .mportant,  at  about 
2000°K,  was  still  increasing  sharply  and  had  t  value  of  greater  than 
2  cal  g-1  deg-1  whereas,  as  Jacobs15  pointed  out,  the  classical  maximum 
for  an  ideal  gas  is  3nR/M  =0,7  where  n  is  the  number  of  atoms,  R  is 
the  gas  constant  and  M  is  the  molecular  weight. 

At  86  kbar,  Campbell,  Davis,  and  Travis1”  have  calculated  a  shock 

temperature  of  1140°K  from  T  =  300  +  Ae/C^,  where  Ae  is  the  increase 

in  internal  energy  across  the  shock.  The  apparent  agreement  with  the 

temperature  calculated  by  us  using  the  Walsh-Christian  method  (Fig.  11) 

is  fortuitous  as  Campbell,  Davis,  and  Travis  used  the  value  of  for 

C  .  If  the  value  of  C  is  used,  the  shock  temperature  calculated  from 
v  v 

this  equation  of  state  becomes  1450°K. 

Also  at  86  kbar,  Mader17  calculated  a  shock  temperature  of  1168°K 

using  the  Walsh-Christian  method,  in  apparent  agreement  with  the  results 

of  similar  calculations  in  Fig,  11.  However  this  agreement  is  also 

fortuitous.  Mader  used  C  instead  of  C  but  the  expected  difference 

p  v 

in  calculated  shock  temperature  was  exactly  offset  by  a  2.5-fold 
difference  in  the  values  used  for  Op/3T)v. 

The  lower  temperatures  calculated  for  nitromethane  using  Cv(T) 
imply  that  the  rate  of  decomposition  of  nitromethane  at  high  pressures 
is  considerably  faster  than  that  expected  from  homogeneous  thermal 
explosion  theory  using  the  rates  measured  at  low  pressures. 

3.  Water 

In  many  ways  water  is  an  anomalous  liquid.  For  this  reason  we 
have  not  attempted  to  estimate  the  variation  of  Cv  with  temperature. 


40 


NITROMETHANE 

TEMPERATURE  —  °K 


PRESSURE  —  kbar 

TA-4101 


FIGURE  11  SHOCK  TEMPERATURE  FOR  NITROMETHANE.  The  input  date  are 
In  Table  2. 
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In  Fig.  12  we  report  the  results  of  calculations  using  the  Walsh- 
Christian  method  and  constant  Cy  and  Op/dTJ^.  The  results  are 
compared  with  those  calculated  by  Rice  and  Walshs  who  assumed  p  to 
be  constant  and  OH/gV)  to  be  a  function  of  pressure  only.  From  the 
analysis  of  the  sensitivity  of  the  Wnlsh-Christian  method  to  and 
Op/dT)  we  have  concluded  that  the  shock  temperature  will  be  very 

sensitive  to  the  value  chosen  for  C  .  The  observed  difference  between 

v 

the  present  results  and  those  obtained  by  Rice  and  Walsh  is  therefore 
regarded  as  lot  significant. 


As  found  for  carbon  tetrachloride,  shock  temperature  is  not 
sensitive  to  the  form  of  the  p-v  Hugoniot ,  although  the  p-v  Hugoniots 
may  be  significantly  different  (Fig.  13). 


G.  Conclusions 

The  present  method  for  calculating  shock  temperature  takes  better 
account  of  the  properties  of  liquids  and  the  greater  dependence  of  shock 
t  mperature  on  than  on  (gp/dT)^.  It  is  therefore  considered  to  be 
an  improvement  on  the  Walsh-Christian  method  and  will  yield  more  realistic 
values  of  shock  temperature  in  liquid  explosives.  This  conclusion 
is  substantiated  by  the  improved  agreement  between  calculated  and 
experimental  temperatures  for  carbon  tetrachloride. 


H.  Future  Work 

Work  on  shock  initiation  is  nearing  completion.  The  new  equation 
of  state  discussed  in  this  report  will  be  used  to  calculate  shock 
temperatures  in  the  dif luoramino  compounds  for  theoretical  calculation 
of  reaction  times. 

In  the  field  of  detonation  failure,  we  plan  to  investigate  the 
stability-structure  problem  in  the  difluoramino  liquids  giving  particular 
emphasis  to  failure-wave  pressure  measurements  and  studies  of  detonation 
wave  structure. 


42 


SHOCK  TEMPERATURE  OF  WATER 


TEMPERATURE  —  °K 


IGURE  12  SHOCK  TEMPERATURE  FOR  WATER.  Comparison  of  the  constant  Cv  model 
with  the  Rice-Walsh5  constant  Cp  model.  The  shock  temperatures  calculated 
using  the  constant  Cv  model  are  not  sensitive  to  the  form  of  the  Hugoniot. 
The  input  data  are  in  Table  2.  See  also  Figure  13. 
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p— V  HUGONIOTS  FOR  WATER 

PRESSURE  —  kbar 


IOC 


DC 

VOLUME 


fl.T 

—  cc  g"1 


FIGURE  13  COMPARISON  OF  SHOCK  HUGONIOTS  FOR  WATER  CALCULATED  FROM  THE 
"UNIVERSAL  HUGONIOTS"  WITH  THE  EXPERIMENTAL  VALUES  MEASURED 
BY  RICE  AND  WALSH5  .  The  input  data  are  in  Table  2.  See  also  Figure  12. 
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Table  6 


COMPARISON  OF  SHOCK  TEMPERATURES  FOR  CARBON  TETRACHLORIDE 
CALCULATED  FROM  EQ.  (7)  BY  TRAPEZOIDAL  EVALUATION  OF  THE  INTEGRAL 
AND  SHOCK  TEMPERATURES  CALCULATED  BY  A  RUNGE-KUTTA  INTEGRATION  OF  EQ.  (6) 


p 

V 

T 

Equation  (7) 
Trapezoidal 

T 

Equation  (6) 
Runge-Kutta 

kbar 

cc  g-1 

°K 

°K 

0 

.631 

298 

298 

29 

.431 

661 

662 

49 

.401 

980 

980 

73 

.381 

1401 

1402 

113 

.361 

2198 

2198 

144 

.351 

2866 

2867 

188 

.341 

3857 

3857 

253 

.331 

5385 

5386 
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Table  7 


INPUT  DATA  FOR  CALCULATING  SHOCK  1 


Constant 

■  cv  fit* 

Mol 

Op/dT) 

Sound 

Specific 

CB/T* 

CC/T 

CD 

CE  X  T 

CF  X  T2 

at  p  =  0 

Units 

_  — 1  ,  -1 
mole  deg  - 

_ _ 

-7 

10  dynes 

10  5  cm 

-1 

cc  g 

-2  .  -l 

-1 

cm  deg 

sec 

Abbreviation  for  Runge-Kutta 

Computer  Program 

CA 

CB 

cc 

CD 

CE 

CF 

CM 

DPDT 

CO 

VO 

Carbon  Tetrachloride 

(1) 

21.7 

0 

0 

0 

0 

0 

153.84 

1.14 

.926 

.631 

(2) 

19.53 

0 

0 

0 

0 

0 

153.84 

1.14 

.926 

.631 

(3) 

21.7 

0 

0 

0 

0 

0 

153.84 

1.254 

.926 

.631 

(4) 

21.7 

0 

0 

0 

0 

0 

153.84 

1.14 

.926 

.631 

(5) 

21.7 

0 

0 

0 

0 

0 

153.84 

1.14 

.926 

.631 

(6) 

21.7 

0 

0 

0 

0 

0 

153.84 

1.14 

.926 

.631 

(7) 

21.7 

-75415.8 

-2109.31 

8.10247 

-8.64548 

1.12516 

153.84 

1.14 

.926 

.631 

*0* 

l 

o 

rH 

X 

x  io~7 

Nltromethane 

(8) 

17,8 

0 

0 

0 

0 

0 

61 

1.637 

1.30 

.884 

(9) 

17.8 

1.23375 

-9956.86 

17.3573 

8.09421  2, 

.24624 

61 

1.637 

1.30 

.884 

X  106 

X  io"3  X 

io-6 

Water 

(10) 

14.07 

0 

0 

0 

0 

0 

18.02 

4.04 

1.48 

1,002 

(11) 

14.07 

0 

0 

0 

0 

0 

18.02 

4.04 

1.48 

1.002 

*  Cv  >  CA  +  (CB/T2) 

+  (CC/T) 

+■  CD  + 

(CE  X  T)  + 

(CF  X  T2) 

/CM 

Source,  of  Data:  Carbon  Tetrachloride :  Cv  (Ref.  9),  cv  fit  (Ref.  18);  (dp/ST)v  (Ref.  9);  Sound  Speed  (Ref.  19); 

Specific  Volume  (Ref.  9);  Hugonlot  1.2,  1.7  (Ref.  2);  1,31,  1,61  (Ref.  8);  1.25,  1.7  (Ref.  this  work, 
arbitrary  variation  of  uy) ;  1.199,  1,672  (Ref.  7) 

N1  tromethane :  Cv  (Ref.  20);  Cv  fit  (Ref.  21);  (dp/dT)v  calculated  from  (dp/dT)v  =  acvc2/v0C  where 
a  la  the  coefficient  of  expansion  and  Cp  is  the  specific  heat  and  constant  pressure  (Re?.  20^; 

Sound  Speed  (Ref.  20);  Specific  Volume  (Ref.  20);  Hugonlot  (Ref.  2) 

Water:  Cv,  (dp/dT)v,  sound  speed,  and  specific  volume  (Ref.  22);  Hugonlot  1.2,  1.7  (Ref.  2);  1.28,  1.58  (Ref.  S 
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[Table  7 


■FULATING  shock  temperatures 


Specific 

Constants  for 

Temp 

Second 

Volume 

First 

Volume 

Temp  i 
Firi 

Volume 

Universal 

at  p» 0 

Point 

Point 

Volui 

at  p  a  0 

Hugoniot 

and  v  «  0 

on 

on 

on 

-1 

U  =  u  e  +  u  u 
s  10  2  p 

o„ 

Hugoniot 

-1 

Hugoniot 

-1 

Hugon: 

o„ 

cc  g 

K 

cc  g 

cc  g 

K 

VO 

U1 

U2 

TO 

XE 

C 

Y[l] 

.631 

1.2 

1.7 

298 

.621 

.631 

298 

.631 

1.2 

1.7 

298 

.621 

.631 

298 

.631 

1.2 

1.7 

298 

.621 

.631 

298 

.631 

1.31 

1,61 

298 

.621 

.631 

298  | 

.631 

1.25 

1.7 

298 

.621 

.631 

298  \ 

.631 

1.199 

1.672 

298 

.621 

.631 

298  j 

.631 

1.2 

1.7 

298 

.621 

.631 

298’ 

.884 

1.2 

1.7 

298 

.874 

.884 

298 

.884 

1.2 

1.7 

298 

.874 

.884 

298 

1,002 

1.2 

1.7 

293 

.809 

.819 

323 

1.002 

1.28 

1.58 

293 

.809 

.819 

323 

Constant  c  In  Flf.  B;  original  In  Fig.  9;  Tables  6,  8,  and  9 

10%  decrease  In  Cy  In  Fig.  9;  Tsbla  B 

10%  Increase  in  (dp/dT)v  in  Fig  9;  Table  8 

Effect  of  p-v  Hugoniot 
Fig.  10  and  Table  9 

Cy(T)  in  Fig.  8 


Constant  Cv  in  Fig.  11 
CV(T)  In  Fig.  11 


Figs.  12  and  13 
Figs.  12  and  13 


work, 

re 

.28,  1,38  (Ref.  23) 
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Table  8 


SENSITIVITY  OP  THE  SHOCK  TEMPERATURE  OP  CARBON  TETRACHLORIDE 


TO  THE  VALUES  OF  Cy  and  (ap/3T)y 


bTs<V0-Vx) 

deg. 

bI/2Cr 

deg. 

■ 

p 

kbar 

V 

cc  g"1 

| 

i 

AT(dp/3T)v 

calc. 

analytically 

atOp/st)v 

obi. 

empirically 

0 

.631 

0 

0 

0 

- 

- 

_ 

29 

.431 

169 

98 

223 

3.6 

1.8 

2.4 

49 

.401 

207 

150 

515 

4.7 

2.4 

3.4 

73 

.381 

233 

208 

978 

6.2 

3.2 

4.5 

113 

.361 

262 

274 

1696 

8.0 

4.1 

5.7 

144 

.351 

277 

300 

2354 

9.8 

5.1 

7.0 

188 

.341 

292 

328 

3335 

12.4 

6.4 

8.3 

253 

.331 

308 

347 

4853 

16.3 

8.4 

9.7 

Input  data  used:  See  Table  7 
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Table  9 


P-v  HUGONIOTS  FOR  CARBON  TETRACHLORIDE  CALCULATED  FROM 
Us  =  +  UjUp 


e  c 

II  II 

1.2 

1.7 

1.25 

1.7 

Hi 

1.199 

1.672 

V 

P 

P 

p 

P 

cc  g"  1 

kbar 

kbar 

kbar 

kbar 

.631 

0 

0 

0 

0 

.431 

29 

32 

31 

28 

.401 

49 

54 

50 

47 

.381 

73 

79 

70 

68 

.361 

113 

122 

103 

103 

.351 

144 

156 

127 

130 

.341 

188 

204 

158 

.331 

253 

275 

202 

In  all  cases  c  =  .926  x  105  cm  sec"1, 
o 
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IV  KINETICS  AND  MECHANISMS  OF  THERMAL  DECOMPOSITION 


(David  S.  Ross  and  M.  E.  Hill) 


A.  Introduction 

Our  last  report1  presented  a  theoretical  model  for  detonations  in 
liquids  in  terms  of  the  types  of  the  chemical  reactions  available  to 
the  material  in  question.  For  NF  materials,  LVD  was  described  as 
resulting  from  an  exothermic  chain  reaction.  HVD  was  proposed  to  be 
due  to  an  exothermic  unlmolecular  dehydrof luorination  for  materials 
such  as  1,2-DP  which  are  structurally  capable  of  undergoing  unlmolecular 
loss  of  HF,  and  some  unspecified  reaction  for  materials  such  as  2,2-DP 
which  cannot. 

There  is  more  to  be  done  with  the  model  in  terms  of  refinements 
and  further  application  to  other  explosive  systems.  However,  we  have 
delayed  further  work  in  this  area  to  allow  time  for  some  laboratory 
work  as  well  as  for  preparation  of  manuscripts  on  NF  chemistry  for 
publication.  We  will  present  below  the  final  results  on  the  very  low 
pressure  pyrolysis  (VIPP)  of  NF  compounds,  as  well  as  our  introductory 
experiments  with  nitropropanes. 

B.  VUP  of  NF  Materials 

The  VLPP's  of  t -but yldifluor amine  (TBD)  and  2,2-DB  have  been 
performed,  and  the  complete  set  of  VLPP  data  on  NF  materials  is  presented 
in  Fig.  14.  These  results  confirm  our  earlier  discussion  of  structure 
and  stability8  in  which  we  suggested  that  there  are  two  factors  which 
act  to  destabilize  the  compound 


They  are,  first,  geminate  substitution  (n>l),  which  lowers  the  activation 
energy  toward  decomposition,  and,  second,  B -substitution,  which  raises 
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FIGURE  14  VLPP  DATA  FOR  NF  COMPOUNDS  AND  DTBP 


the  A  factor  for  decomposition.  2,2-EB  is  destabilized  through  both 
processes,  and  thus  as  we  had  predicted,2  is  the  least  thermally  stable 
material  we  have  studied. 


A  discussion  of  the  detailed  mode  of  decomposition  of  vicinal 
bisdifluoramines  entails  consideration  of  the  reaction  from  both  directions. 
Thus  the  fact  that  propylene,  butene-2,  and  Isobutylene  are  produced 
as  the  overwhelmingly  major  products  in  the  VLPP  of  1,2-DP,  2,3-DP,  and 
1BA  respectively  suggests  that  we  are  observing  the  reverse  of  the 
addition  reaction  of  NaF4  to  these  olefins.  This  addition  has  been 
studied  by  Dijkestra,  Kerr,  and  Trotman-Dickensonsa>b>c  (DKT) ,  and  a 
detailed  mechanistic  path  for  the  reaction  sequence  based  on  both  that 
work  and  ours  is : 


1 

F2N-C-C-NFa  +  M  i*  FaN-C-C-NFa*  +  M 
-1 

*  2 

FaN-C-C-NFa  «*  FjN-C-C"  +  -NFa 
-2 
3 

FaN-C-C-  +  M  i*  FaN-C-C.  +  M 
-3 

*  4 

FaN-C-C-  FaN-  +  C=C 
-4 

5 

M  +  2.NFa  ;*  N?F4  +  M 
-5 


(1) 

(2) 

(3) 

(4) 

(5) 


in  which  *  denotes  a  vibrationally  excited  state  and  M  is  a  collision 
partner.  Each  step  is  written  as  reversible,  however,  depending  upon 
the  method  and  direction  of  study,  certain  steps  proceed  only  in  one 
direction  in  a  practical  sense.  Thus,  for  example,  Eq.  (4)  in  the 
case  of  the  work  by  DKT  is  truly  reversible,  whereas  for  this  work, 
step  4  is  very  rapid  but  step  -4,  because  of  the  low  pressure,  proceeds 
essentially  at  zero  rate. 

Figure  15  provides  a  schematized  reaction  coordinate  for  the 
sequence,  and  embodies  the  fact  provided  by  DKTsb  that  the  activation 
energies  for  3  and  -4  are  the  same  to  within  2  kcal/mole,  and  thus  the 
C-NFa  bond  strength  in  FaN-C-C-  is  small.  They  fall  in  the  range 
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8-16  kcal/mole  for  the  olefins  studied,  and  the  value  used  in  Fig.  15 
of  ~14  kcal/mole  is  for  propylene.  Since  kj  is  small,  v,  (the 
velocity  of  step  2)  »  v_!,  and  thus  step  1  is  rate  controlling.  The 
experimental  activation  energy  for  vicinal  blsdif luoramlne  (~  55  kcal/mole) 
applies  therefore  to  step  1.  The  heat  value  of  ~  20  kcal/mole  for  -5 
is  from  the  work  of  Piette,  Johnson,  Booman,  and  Colburn.4 


FIGURE  It).  REACTION  COORDINATE  DIAGRAM  FOR 
FOR  THE  REACTION: 

C-C  +  n2f45=  f2n  -  C  -  C  -  nf2 


DKT  showed  that  in  their  direction  of  study  the  slow  step  is  -4, 
and  while  Fig.  15  suggests  that  both  -4  and  -5  are  slow,  in  fact 
v_5  »  v_4.  The  reasoning  here  points  to  the  need  to  be  careful  in 
interpreting  reaction  coordinate  diagrams  in  that  they  show  the  enthalpy 
but  not  the  entropy  factors  involved  in  the  process.  Thus  at  the 
temperatures  and  pressures  of  the  DKT  work 

i 

k_5  =  10‘ 3- 5  exp  (-15, 200/RT)  l/(mol  sec) 


under  DKT  conditions  of  a  mean  temperature  of  90°C  and  pressures 
in  the  range  20  to  600  torr,  is  in  its  second  order  region,5  and  this 
expression  is  derived  from  k  =  1.5  x  I0la  T  0>  5  exp  (-15, 200/RT)  liters/ 
(mol  sec)5  for  M  =  argon  and  with  T  =  400°K. 
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while 

k_4  =  10*  •«  exp  (-13 i 700/RT)  1/(bo1  sec) 

The  activation  energy  terma  are  comparable,  but  the  preexponential 
terma  differ  by  about  aix  powera  of  ten.  In  fact,  under  the  DKT  experi¬ 
mental  conditiona,  Eq.  (5)  la  a  true  equilibrium. 

It  la  of  intereat  to  conaider  the  eventa  along  the  plateau  of 
Eq.  (2)  in  Fig.  IS.  Here  the  vibratlonally  excited  apeciea  F,N-C-C-NFa* 
is  undergoing  rapid  Internal  exchanges  of  energy  among  the  varloua 
modes  available,  and  while  in  most  cases  a  C-N  bond  is  ultimately 
broken,  there  is  another  option  available  to  the  molecule.  It  has 
recently  been  shown  that  vibratlonally  excited  CH,NFa  loses  HF  yielding 
HCN."  Thus,  it  is  conceivable  that  such  a  process  la  taking  place 
here  in  the  case  of  1,2-DP,  for  example, 

CHa-CH-CHaNFa*  -I-J-  3HF  +  CHaCCN 
NFa  (F)NF 

which  would  explain  the  presence  of  N-fluorlminoproplonitrile  in  the 
product  mixture  from  that  compound.  The  SiF«  would  then  arise  from 
the  well  known  reaction  of  HF  with  glass 

4HF  +  SiF,  -*  2Ha0  +  SiF4 

Lesser  products  from  the  VLPP  of  1,2-DP,  such  as  CHaCNa,  rcHaCN, 
and  CHjNF,  car  be  explained  by  considering  the  fact  that  the  effects 
of  pressure  on  the  DKT  kinetics  of  HaF4  addition  to  olefins  suggested 
that  the  vibratlonally  excited  species  FNa-C-C- *  has  some  finite  lifetime 
(as  opposed  to  being  a  transient  intermediate).  In  view  of  our 
isolating  a  material  with  a  C-F  bond,  we  propose  that  there  is  a  route 
available  to  the  oxcited  radical  in  addition  to  -3  and  4,  namely 
internal  migration  of  fluorine  from  nitrogen  to  carbon. 
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JM 

-3 


FaN-C-C 


FjN-  +  C=C 

FN-C-CF 


This  process  has  adequate  driving  force  as  it  is  some  35-45  kcal/mole 
exothermic,  and  since  DKT  operated  at  pressures  in  the  20  to  600  torr 
range  while  this  work  was  performed  at  less  than  10" 3  torr, 

DKT  This  work 

V.J  »  V_3 


Thus  6  is  observable  in  our  work  whereas  it  was  not  in  that  by  DKT. 

Thus,  for  the  VLPP  of  1,2-DP  the  following  scheme  is  proposed  to 
account  for  the  minor  products  outside  of  the  already  explained 
N-  f  luoriminoproplonit  rile : 


CHjCH-CH.* 

I 

NF, 

~F 


CH3CN  FCHjCN 

The  methyl  radicals  are  accounted  for  in  the  observation  of  CH3NFa 
and  HCN6 
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2HF  +  HCN 


CHj  +  NFa 


CHjNFj 


The  same  process  for  the 


fragment  CH3CHCHaNFa* 


(7) 


CH,CHCH,NF„* - Z?  >  CH-CHFCH- 

'  I 

NF 


.1* 


CH3CHF  +  HaC=NF 


-HCN 


-HF 


(8) 


must  take  place,  however,  the  product  HaC=NF  has  not  been  detected, 
and  the  small  amount  of  HCN  observed  can  be  accounted  for  by  Eq,  (7). 
One  explanation  for  this  observation  is  that  the  secondary  radical 
in  8  is  2-3  kcal/mole  more  stable  than  the  primary  radical  in  6,  7 

and  thus  the  latter  should  undergo  rearrangement  more  rapidly  than  the 
former. 


Of  course,  it  is  possible  to  explain  the  formation  of  all  of  the 
minor  products  by  wall  reactions.  In  fact,  we  have  evidence  that  1,2-DP 
does  dehydrofluorinate  readily  through  a  surface-catalyzed  process  to 
yield  N-fluoriminopropionitrlle.  *  However,  since  there  are  plausible 
homogeneous  routes  leading  to  the  observed  products,  we  prefer  to  suggest 
them  as  the  route  of  formation. 


C.  VLPP  of  Nltro  Materials 

VLPP'a  of  both  1-  and  2-nit ropropane  (1-NP  and  2-NP)  have  been 
reported  by  Spokes  and  Benson.  B  We  are  repeating  the  work  for  the 
sake  of  completion  as  a  prelude  to  our  study  of  the  dinitropropanes. 
Spokes  and  Benson  using  a  9,000  collision  vessel  found  that  the  nltro- 
propanes  decompose  with  the  following  high  pressure  rate  parameters 

kl-NP  =  10*,-5±0-«~<«*±*>/e  sec-1 
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and 


_  jqII.  *±o.  a-(4Qfco.5)/0 


k 

2-NP 


sec” 1 


(9  =  2.303  RT  x  10J  kcal/mole)  via  a  S-center  transition  state, 
eliminating  HNO,  in  a  concerted  fashion  and  yielding  propylene.  Thus, 
for  example,  the  VLPP  of  2-NP  proceeds  as  follows: 


CHjCHCH,  - 
NO, 


H,C;£'CH’ 


* 


CH3CH=CH, 


+  HONO 


We  have  studied  1-NP  in  a  100  collision  vessel  and  have  found  a 
Ti/a  °*  665°C.  Propylene,  as  expected,  is  the  single  observed  product. 
Wh.1  le  HNO,  cannot  be  seen  under  our  gas  chromatographic  conditions, 
it  is  most  likely  present,  2-NP  is  currently  being  studied. 

The  study  of  the  dinit ropropanes  (DNP)  is  expected  to  be  complicated 
by  the  fact  that  these  materials  have  rather  low  vapor  pressures  at  room 
temperature.  Thus  at  25°c  1,1-DNP  and  2,2-DNP  have  vapor  pressures 
well  below  1  torr,  and  a  vapor  pressure  convenient  for  VLPP  operation 
of  5  to  10  torr  is  attainable  only  at  temperatures  of  70  to  80°C.,D 
The  sample  volume  of  the  reactor  will  have  to  be  heated  to  this  tem¬ 
perature  to  allow  the  study  to  proceed,  and  we  are  currently  considering 
various  possible  ways  of  doing  this. 

The  thermal  decomposition  of  2,2-DNP  from  175  to  210°C  has  been 
studied  by  Flournoy11  and  in  that  work  he  observed  equimolar  proportions 
of  acetone,  NO,  and  NO,  as  products.  The  proposed  mechanism  involved 
initial  scission  of  the  CN  bond 


(CH,),C(N0,),  -  NO, •  +  (CH3),CN0,  (1) 

followed  by  formation  of  an  Intermediate  nitronitrite  by  either 
recombination 

.  0110 
(CH,),CN0,  +  NO,’  -  (CHj),C  (2) 

Nf0, 
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or  by  displacement  on  another  molecule  of  2 , 2-DNP 

.  •  yNO, 

(CH3)aCNOa  +  <CHa)aC(NOa)a  ->  (CH,),C  +  (CH,)tCNO.  (3) 

OHO 

The  final  step  In  the  formation  of  acetone  la  proposed  as  the  rapid 
decomposition  of  the  nitronitrite 


(CH3)aC  1  -  (CHj)jC=0  +  NO-  +  NO 
0N0 

The  rate  of  disappearance  of  2, 2-DNP  followed  the  rate  law 


(4) 


k 


2,  2-DNP 


10ia.s-so/9  see-i 


which  was  applied  to  the  slow  step  (1) 


A  comparison  of  the  constants  and  k>|M|p  shows  that  for 

the  case  of  geminate  materials  C-N  scission  will  be  faster  than  HONO 
elimination  above  300°K.  The  VLPP  study  of  2,2-DNP  should  prove 
Interesting  since  VLPP  operation  is  at  such  low  pressure,  the  rate  of 
either  step  2  or  step  3  should  be  vanishingly  small,  and  thus  we  should 
see  the  products  resulting  from  the  decomposition  of  the  radical 
intermediate,  (CH3)2CNOa  .  For  this  process  there  are  two  possible 
pathways.  First,  by  analogy  to  the  decomposition  of  t-butyl  radical 
in  our  reactor, 12  the  intermediate  could  lose  a  hydrogen  atom  to 
yield  2-nitropropene 

CH, 

I  * 

•C-NO,  -  H-  +  CH3C=CHa 
CHS  N0a 


Second,  an  intramolecular  oxygen  rearrangement  may  take  place  to  yield, 
subsequently,  acetone  and  NO. 

CH,  oe  CH j 

•c-f'  9-c-f;  -  (CH3)2C=0  +  NO 

<*,  °  <1h3  0 
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This  second  possibility,  by  the  way,  coming  after  initial  loss  of 
N0a,  would  provide  equimolar  quantities  of  acetone,  NO,  and  NOa,  and 
thus  is  an  alternate  explanation  of  Flournoy's  results. 

D.  Future  Work 

Woolfolk  has  presented  an  approach  to  correlating  LVD  gap 
sensitivity  to  thermochemical  stability  (see  Section  II).  His  reaction 
heat  values  (AHr)  are  defined  by  the  process  whereby  the  liquid 
explosive  goes  to  thermodynamically  favored  end  products  such  as  HF, 

C,  and  Na .  While  at  this  point  such  a  definition  for  AHr  is  all  that 
is  reasonably  available,  it  is  not  clear  that  in  an  LVD  process  such 
products  are  indeed  formed.  For  example,  since  the  temperature  in 
LVD  may  never  exceed  1000°C,  it  is  possible  that  less  stable  (in  a 
thermodynamic  sense)  materials  such  as  olefins,  imines,  and  nitriles, 
are  "trapped,"  in  that  kinetically  they  cannot  proceed  to  the 
aforementioned  end  products. 

Were  some  measure  of  the  degree  of  the  thermodynamic  inefficiency 
known  for  the  LVD  process,  a  correction  could  be  applied  to  the  AH^ 
values,  and  a  firmer  correlation  would  be  the  result.  In  this  regard 
we  have  discussed  the  possibility  of  carrying  out  LVD  on  a  small  scale 
in  some  sort  of  enclosure  such  that  the  products  of  the  explosion 
could  be  trapped  and  characterized.  Particularly  interesting  would 
be  the  possibility  of  finding  products  in  these  proposed  experiments 
that  have  been  also  seen  in  the  VUP  studies  on  the  same  compounds. 
Thus,  for  example,  the  presence  of  FCHSCN  as  one  of  the  end  products 
of  the  LVD  of  1,2-DP  would  verify  our  supposition1 3  that  the  highly 
exothermic  migration  of  fluorine  from  nitrogen  to  carbon  as  observed 
by  VLPP  plays  an  important  role  in  the  detonation  process. 
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and  geminate  difluoramino  groups  on  the  ease  of  initiation  and  on  the 
environmental  stability  of  compounds  containing  them.  A  similar 
approach  has  begun  with  isomeric  nitro  compounds. 

Low-Velocity  Detonation 

The  shock  wave  Interactions  that  can  lead  to  the  initiation  of 
low-velocity  detonations  (LVD)  are  examined  with  particular  emphasis 
on  cavitation.  The  two  forms  of  cavitation,  gas  inclusion  and  pure 
vaporization,  are  discussed.  The  latter's  relationship  to  the  shock- 
wave  interactions  in  the  LVD  gap  sensitivity  teats  is  examined.  We 
believe  that  pure  vaporization  cavitation,  which  cannot  be  removed 
by  degassing,  is  present  under  gap  test  conditions.  Experimental 
studies  in  the  gap  test  arrangement  show  the  presence  of  cavitation  in 
qualitative  agreement  with  theoretical  prediction. 

Tile  universal  Hugoniot  for  liquids  has  been  reevaluated  using  a 
temperature  correction  for  the  sonic  velocity,  c0.  The  new  equation  is 
improved  slightly  over  the  one  formerly  used.  This  concluded  the 
effort  in  this  area, 

A  trace  dye  method  for  following  the  flow  field  behind  a  shock 
wave  was  studied.  We  concluded  that  this  method  was  not  suitable  at 
present  for  our  studies  and  work  on  it  will  be  discontinued. 

Experiments  on  the  effect  of  confinement  wall  thickness  on  the 
LVD  gap  sensitivity  are  discussed  and  preliminary  charge  diameter  studies 
are  reported.  A  summary  of  the  effects  of  confinement  on  the  LVD 
gap  sensitivity  is  presented. 

Some  of  the  physical  and  chemical  properties  of  liquid  explosives 
are  correlated  with  the  LVD  gap  sensitivity.  This  correlation  is 
remarkably  good  considering  its  simplicity.  The  evidence  that  LVD 
initiation  may  be  a  gas  phase  phenomenon  is  discussed.  This  further 
supports  the  contention  that  LVD  may  be  related  to  the  formation  of 
cavities  in  the  liquid  explosive. 

Physics  and  Chemistry  of  Detonation 

A  new  method  for  calculating  shock  temperatures  in  liquid  explosives 
has  been  developed  from  the  Walsh-Christian  method  for  metals.  The 
Walsh-Christian  method  assumes  that  Cv  and  Op/^T)v  are  constant.  The 
basis  of  the  present  method  is  a  variational  analysis  which  shows  that 
shock  temperature  is  more  sensitive  to  variations  in  Cv  than  to  variations 
in  Op/5T)v,  Thus  Cy  is  assumed  to  be  a  function  of  temperature  and 
0P/9T)V  is  assumed  to  be  constant.  It  is  considered  likely  that  Cv 
will  Increase  with  temperature  due  to  an  increase  in  the  vibrational 
heat  capacity.  The  simplest  assumption  is  that  C  of  the  liquid  will 
increase  as  Cv  of  the  ideal  gas.  This  is  the  basis  of  the  CV(T)  model. 

Shock  temperatures  of  carbon  tetrachloride  and  nitromethane  have 
been  calculated  using  both  constant  Cy  and  C  (T).  Shock  temperatures 
for  carbon  tetrachloride  calculated  using  the  CV(T)  model  are  in  better 
agreement  with  experimental  brightness  temperatures  than  those  calculated 
using  the  Walsh-Christian  constant  Cv  model.  Shock  temperatures  for 
nitromethane  are  compared  with  those  of  other  workers. 
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Shock  temperatures  in  water  were  calculated  using  constant  Cv  and 
OP/£T)v.  The  temperatures  are  about  15 %  less  than  those  calculated 
by  Rice  and  Walsh  using  a  model  based  on  constant  Cp.  This  difference 
is  not  regarded  as  significant  compared  with  the  uncertainties  in  the 
variation  of  Cv  and  Cp  with  temperature  and  pressure. 

Kinetics  and  Mechanisms  of  Thermal  Decomposl'tlon 

The  VU*P  studies  of  all  the  NF  compounds  we  have  had  at  hand  have 
been  completedi  and  the  accumulated  data  suggest  a  model  of  unimolecular 
decomposition  in  which  the  NF2  compound 

C.  -  C  (NF.) 

@  a  a  n 

is  destabilized  by  both  a~NFa  substitution  due  to  a  decrease  in  Ea>  and  by 
0-alkyl  orp-NF2  substitution  due  to  an  increase  in  the  A  factor.  Since 
both  factors  are  operative  in  the  case  of  2t2~DB>  it  is  thermally  the 
least  stable  NF  material  we  have  studied. 

A  detailed  reaction  coordinate  diagram  may  be  constructed  for  the 
system  F2N-C-C-NF2  £  C=C  +  N2F4  on  the  basis  of  both  our  work  and  that 
by  Trotman-Dlckenson.  It  is  possible,  on  the  basis  of  the  diagram,  to 
explain  the  side-products  observed  in  the  VLPP  of  1,2-DP. 

The  VLPP  of  1-nit ropropane  has  been  studied.  The  products  are 
propylene  and  presumably  HONO,  and  the  T, /2  value  is  665°C. 
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